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PREFACE 


This  is  the  third  in  a  continuing  series  of  Memoranda 
on  the  development  of  a  mathematical  simulation  of  human 
blood  biochemistry.  The  first--Dantzig,  et  al.  [lT--was 
a  feasibility  study  of  the  methods  for  mathematical  anal** 
ysis  and  computer  simulation;  the  second- -DeHaven  and 
DeLand  [2]--added  some  details  of  hemoglobin  chemistry 
and  carefully  described  the  behavior  of  that  system  under 
various  conditions. 

The  subsequent  body  of  critical  comment  and  addi¬ 
tional  laboratory  data  enabled  further  elaboration  and 
definition  of  the  model.  Professor  F.J.W.  Roughton,  F.R.S., 
Cambridge  University,  particularly  improved  the  hemoglobin 
chemistry.  In  addition,  the  quantitative  aspects  of  Lae 
plasma  proteins  have  been  reviewed  and  their  buffering 
properties  incorporated. 

The  present  model  therefore  approaches  more  closel> 
the  complexity  of  the  real  blood  system.  Properties  of 
the  mathematical  model  such  as  gas  exchange,  buffering, 
and  response  to  chemical  stress  in  the  steady  state  are 
practically  indistinguishable  from  those  properties  of 
real  blood  within  the  limits  of  our  current  validation 


program. 


r 


The  program  consists  of  two  parts:  detailed  chemical 
analysis  of  human  blood  under  a  variety  of  chemical  stresses; 
and  further  elaboration  of  the  model,  e.g. ,  by  including 
the  phosphate  system,  double-valent  Ion  binding  by  proteins, 
and  red  cell  organic  constituents  ocher  than  hemoglobin. 
Laboratory  experiments  validating  the  model  will  be  pub¬ 
lished  in  the  doctoral  thesis  of  Eugene  Magnler,  M.D. , 

Temple  University  Medical  School. 

The  research  done  with  Dr.  F.J.W.  Roughton  was  sup¬ 
ported  by  the  Department  of  Health,  Education,  and  Welfare 
Grant  HE-08665  (National  Institutes  of  Health). 
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SUMMARY 


This  Memorandum,  one  of  a  series  on  the  modeling  of 
blood  biochemistry,  describes,  and  examines  the  consequences 
of,  the  classical  analysis  of  human  blood.  Classically,  the 
macroscopic  features  of  healthy  blood--the  principal  £1  i id 
and  electrolyte  distributions--can  be  said  to  be  in  a 
hemostatic  steady  state  under  the  influence  of  several 
interacting  constraints,  e.g.,  the  osmotic  effects  of  the 
fixed  proteins,  the  (fixed)  charge  of  the  proteins  and 
'the  neutral  electrostatic  charge  condition,  and  the  active 
cation  pumps.  Chapter  I  examines  these  constraints  from  the 
vantage  of  a  theoretical  model,  and  shows  them  to  be  suf¬ 
ficient  explanation  for  the  steady  state;  whether  they  are 
also  necessary  has  not  been  demonstrated. 

Chapter  II  considers  the  incorporation  of  the  micro¬ 
scopic  properties  of  the  proteins  in  the  model,  particularly 
their  buffering  behavior.  This  work  is  ot  complete,  owing 
principally  to  the  lack  of  firm  data  on  hemoglobin. 

Chapter  III  discusses  the  validation  of  the  mathematical 
model  and  the  consequences  of  the  previous  biochemical  struc¬ 
tural  detail,  testing  the  model  under  various  conditions-- 
changes  in  gas  pressure  and  pH,  and  additions  of  a  few 
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chetnicfll  stresses.  It  compares  several  validation  experi¬ 
ments  with  the  literature.  The  Memorandum  shows  that, 
except  for  certain  interesting  discrepancies,  the  model 
satisfactorily  agrees  with  the  literature. 
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Chapter  I 

A  PRELIMINARY  MODEL  OF  BLOOD  BIOCHEMISTRY 


1 .  INTRODUCTION 


This  is  the  third  in  a  continuing  series  of  Memoranda 
detailing  the  development  of  a  mathematical  simulation  of 
human  blood  biochemistry.  The  f irst--Dantzig ,  et  al.  [ 1  ]  — 
was  a  feasibility  study  of  a  method  for  mathematical  formu¬ 
lation  and  computer  solution  of  a  biochemical  simulation 
problem.  An  example  of  a  "simple"  problem  is,  given  tbe 
basic  molecular  components  (reactant  species)  of  a  chemical 
milieu  in  a  single  phase,  compute  the  concentration  of  all 
species  in  the  equilibrium  mixture  for  which  the  chemical 
reaction  coefficients  and  mass  action  constants  are  also 
given . 

Such  a  system  may,  however,  become  complex  if  many  re¬ 
actions  are  interrelated,  if  some  of  the  molecules  in  the 
milieu  (such  as  proteins)  are  themselves  complex,  or  if 
the  system  is  multiphasic,  multicompartmented ,  or  alive. 

These  complexities  complicate  the  comput  *t.ion;  in  particular, 
for  viable  systems,  one  speaks  of  the  computations  for  a 
supposed  "steady  state"  of  the  living  organism,  rather  than 
the  "equilibrium  state." 
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The  Dantzig  Memorandum  introduced  the  computer  method 
for  minimizing  the  Gibbs  free  energy  function  to  compute 
the  (unique)  steady  state  of  a  multicompartmented  blood 
system.  The  procedures  are  described  in  Ref.  1,  and  de¬ 
veloped  further  in  Shapiro  and  Shapley  [3],  Shapiro  [4], 
and  Clasen  [5].  DeHaven  and  DeLand  [2]  (the  second  Memo¬ 
randum  on  blood  simulation)  explore  the  computer  method  for 
determining  a  steady  state  in  the  three-compartment  model 
(lung  gas.  plasma,  and  red  cell),  emphasizing  the  properties 
of  the  chemical  systems  themselves.  Even  that  rudimentary 
model  had  many  important  characteristics  of  a  blood  system, 
e.g.,  the  proper  distribution  of  water  and  electrolytes, 
and  a  reasonable  representation  of  normal  physiologic 
functions . 

This  Memorandum  extends  the  previous  results ,  docu¬ 
menting  them  in  more  detail,  and  initiates  the  validation 
of  an  improved  model  of  the  blood  chemistry.  Chapter  I 
analyzes  the  biochemical  structure  of  the  two-compartment 
blood  system.  The  purpose  of  this  chapter  is  to  elucidate, 
with  a  rudimentary  blood  model  such  as  that  used  in 
Dantzig  [1],  the  conventional  roles  of  the  fixed  proteins, 
the  neutral  electrostatic  charge  constraints,  and  the 
active  cation  pumps  in  determining  the  major  physical 


characteristics  of  hemostatic  blood.  While  these  con¬ 
ventional  roles,  as  is  well  known,  yield  a  sufficient 
general  explanation  for  the  fluid  and  electrolyte  distri¬ 
bution  of  normal  blood,  a  mathematical  model,  which  permits 
experiments  impossible  in  the  laboratory,  shows  their 
detailed  interplay. 

Chapter  II  describes  a  special  mathematical  procedure 
for  representing  large  chemical  systems  and  discusses  the 
representation  of  the  essential  proteins.  Finally,  Chap. 
Ill  discusses  a  model  of  the  respiratory  biochemistry  of 
the  blood,  using  the  previously  developed  details,  and 
shows  some  of  the  model’s  initial  validation  procedures. 


ISOLATING  AND  DEFINING  THE  PHYSIOLOGICAL  SYSTEM 


The  blood,  contained  in  the  vascular  system,  may  be 
regarded  as  a  uniform  and  continuous  subsystem  of  the  body. 
Of  course,  it  changes  chemically  in  the  circuit  around 
the  body,  visibly  changing  color.  However,  conceptually 
isolating  the  blood  as  a  subsystem  and  considering  only 
the  biochemistry  of  the  respiratory  function,  this  sub¬ 
system  exhibits  the  following  properties: 

a)  The  subsystem  has  natural  boundaries;  its  extent 
is  well-defined.  It:  is,  for  practical  purposes, 
contained  within  the  vascular  system.  Important 
and  interesting  exceptions  to  this  sometimes 
occur — e.g.,  protein  may  pass  through  the  capil¬ 
lary  walls  to  be  returned  to  the  vascular  system 
by  lymph  flow. 

b)  The  function  of  the  bounding  membranes  can  be 
defined  and  the  substances  crossing  these  mem¬ 
branes  analyzed--i. e. ,  the  inputs  to  and  losses 
from  the  subsystem  are  measurable.  The  actual 
mechanisms  and  biochemical  processes  within  the 
membranes  themselves  are,  for  the  most  part, 
still  hypothetical,  but  only  the  net  results  of 
their  functions  are  necessary  to  simulations 

of  the  steady-state.  Transport  of  substances 
across  these  bounding  membranes  occurs  prin- 
cipally--and  perhaps  only--in  capillary  beds 
either  at  the  lung  or  other  body  tissues. 

c)  Biochemical  analyses  of:  the  blood  can  be  made 
easily  and  quickly  and,  with  care,  thoroughness 
and  accuracy  is  possible.  Furthermore,  samples 
from  this  subsystem  ar,>  available.  Thus,  for 
example,  conservation  of  mass  equations,  funda¬ 
mental  to  the  simulation,  can  be  substantiated. 
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d)  The  subdivisions  of  the  subsystem  on  the  gross 
level  are  clear,  e.g.,  the  red  cells  taken  to¬ 
gether  as  a  compartment  constitute  a  subdivision 
of  whole  blood.  However,  the  exterior  environ¬ 
ment  of  the  vascular  system  must  also  be  examined 
for  a  model,  since  this  environment  largely  de¬ 
termines  the  transfer  of  substances  across  the 
bounding  membranes.  Except  in  the  interior  of 
certain  organs,  the  environment  of  the  vascular 
system  is  interstitial  fluid.  However,  in  the 
lung  (and  perhaps  the  brain)  the  layer  of  inter¬ 
stitial  fluid  is  negligibly  thin.  In  the  lung 

we  may  consider  the  environment  of  the  capillary 
bed  to  be  alveolar  sac  gas. 

e)  Ihe  physiological  function  of  respiration  for 
the  subsystem  "blood"  is  well-defined;  its 
purpose  is  clear.  The  biochemistry  involved  in 
accomplishing  this  purpose  is,  on  the  gross 
level,  relatively  well-understood.  However,  on 
the  level  of  chemical  detail,  many  unsolved 
problems  still  exist.  For  example,  the  explana¬ 
tions  of  the  "chloride  shift"  phenomenon,  or 

the  "Bohr  shift"  phenomenon,  or  the  "sodium  pump" 
phenomena  are,  among  many  other  such  problems, 
still  subjects  of  research. 

The  chemical  composition  of  blood  varies,  of  course, 
according  to  the  point  in  the  body,  the  age  and  sex  of 
the  individual,  his  response  to  disease,  etc.  An  average, 
resting,  young,  adult  male  human  is  the  standard  for  this 
exposition.  The  relevant  data  come  from  man;  sources,  but 
particularly  Dittmer  (Ref.  6,  Blood  and  Other  dy  Fluids) 


and  Spector  (F  f.  7,  Handbook  of  Biological  Data). 
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Table  I  summarizes  these  data  for  three  standard 
types  of  blood  for  the  average,,  young,  adult  male.  Table 
II  shows  the  distribution  of  the  substance  from  Table  I 
for  standard  or  average  artericl  blood  between  plasma  and 
red  cell  compartments.  These  composition  tables  are 
generally  concerned  with  electrolytes,  water,  and  gases. 
This  chapter  does  not  detail  the  amino  acids,  lipids, 
carbohydrates,  organic  phosphates,  or  other  secondary 
constituents.  It  mentions  only  the  constituents  which 
play  a  direct  role  in  respiration. 

The  plasma  proteins  and  impermeable  macromolecules, 
principally  albumin  and  the  various  globulins,  amount 
to  about  8.7  x  10  moles  per  liter  of  blood.  At  pH  7.4, 
the  average  change  per  molecule  for  the  proteins  is  from 
-1.0  to  -20,  as  can  he  determined  by  examining  the  titra¬ 
tion  curves  for  albumin  and  globulin  [3\  For  this  first 
model  where  the  proteins  are  not  yet  buffering,  the  value 
of  negative  equivalents  per  mole  will  be  fixed  at  -10. 
Thus,  the  number  of  negative  equivalents  of  plasma  pro¬ 
teins,  used  to  determine  that  the  solution  is  electrically 
neutral,  is  8.7  x  10  \ 

The  reo  cell  proteins  present  a  much  more  complicated 
situation.  Of  course,  hemoglobin  activity  is  crucial  to 


Table  I 


STANDARD  BLOOD  COMPOSITION  FOR  RESTING  YOUNG  ADULT  MALE 


Item 

Venous  Blood, 
Large  Vein 

Mixed-Venous , 
Pulmonary 
Artery 

Arterial , 
Pulmonary 
Vein 

PH 

7.37 

7 . 3£ 

7.39 

P°2 

65.0  mm 

40.0  mm 

100.0  mm 

pCC^ 

46.0  mm 

46.0  mm 

40.0  mm 

PN2 

573.0  mm 

573.0  am 

573.0  mm 

HCO^ 

22.3  meq/L 

22.0  meq/I 

21.0  meq/L 

70  Sat. 

91.0 

73.0 

97.0 

h20 

46.50  moles/L 

(b) 

(b) 

Cl' 

80.05  meq/L 

Na+ 

84.82 

K+ 

45  05 

Ca^ 

3.210  " 

++ 

Mg 

3.225 

soT 

1.050  " 

4 

HPO" 

4 

2.860  " 

Lact Ic 

2.033  " 

nhT 

4 

0.869  " 

Urea 

3.142  mmoles/L 

Glucose 

3.666  " 

PI  PR'10 

0.870  " 

PC  pr"x 

10.57 

Hb4 

2.2725  " 

A 

See  notes  In  text. 


All  the  following  items  in  this  column  do  not  change 
with  the  type  of  blood,  or  iheir  change  has  not  been  measured, 
or  tht  change  is  not  significant  at  this  s^age  of  the  simulatio 
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1.3.15 

5.26 

75.4 

6.11 


Plasma 

Red  Cell 

7.39 

7.19 

6.33x10  ^  noles 

6.43x10  ^  moles 

7.01xl0‘4  " 

4. 76xl0~4  " 

2.20xl0‘4  " 

2.20xl0"4  " 

28.50 

lb.  0  " 

56.65  meq 

23.40  meq 

76.45  " 

8.37  " 

2.310  " 

42.75  " 

2.86  " 

0.344  " 

0.935  " 

2.295  " 

0.680  ” 

a 

1.88  " 

1.30  M 

— 

0.054  " 

13.836  " 

6.025  M 

2.00  mmoles 

1.73  " 

1.40  mmoles 

0.87  " 

10.568 

2.2725  " 

Indicates  the  value  is  either  not  reported 
or  is  reported  in  a  form  not  pertinent  to  this  representation 
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respiration,  carrying  CO2  as  well  as  oxygen  in  combina¬ 
tion.  In  more  complex  models  (Chap.  Ill),  the  buffering 
activity  of  hemoglobin  also  pertains.  This  preliminary, 
illustrative  model  treats  hemoglobin  simply  as  a  protein 
combining  with  oxygen  and  contributing  to  the  osmolarity 
and,  hence,  hematocrit  of  the  blood. 

However,  the  red  cells  also  contain  miscellaneous 
proteins,  complex  phosphates,  and  other  anionic  species 
(excluding  hemoglobin)  generally  represented  in  the  liter¬ 
ature  [7]  by  X  .  This  symbol  indicates  an  undeterminea 
anion  residue  necessary  to  make  the  red  cell  milieu 
neutrally  charged  and  to  make  the  osmolarity  of  the  red 
cell  milieu  equal  to  that  of  plasma.  Chapter  III  examines 
this  assumption  of  equal  osmolarity,  but  the  preliminary 
model  uses  the  undetermined  anion  residue  to  obtain  the 
proper  hematocrit  and  a  neutral  e'<ctrical  charge. 

Neither  the  distributions  of  glucose  and  urea  nor  of 
the  various  sulphates  and  phosphates  between  plasma  and 
red  cell  have  been  carefully  detailed  for  this  preliminary 
model.  The  diffusion  of  glucose  through  the  red  cell 
membrane  is  presumably  an  intricate  process,  while  urea 
is  thought  to  pass  freely  through  the  membrane.  Table  II 
reflects  neither  of  these  facts.  The  species  represented 


n. 


-10- 


by  SO^  and  HPO^  are  the  acid-aolable  sulphates  and  phos¬ 
phates  reported  in  the  literature.  Table  Ii  does  not 
indicate  that  these  substances  may  also  be  intricately 

bound  to  the  proteins  of  red  cell  and  plasma.  Similarly, 

1 1  1 

Ca  and  Mg  may  be  bound  by  protein.  Chapter  III  will 
discuss  these  conditions. 

Most  of  the  Ha  and  a  larg'  part  of  the  Cl  are  in 
plasma,  while  most  of  the  K+  is  in  the  red  cell  compart¬ 
ment,  a  consequence  of  the  "sodium  pump"  and  the  Gibbs- 
Donnan  phenomenon.  The  distribution  of  the  species  H^O 
in  Table  II  requires  some  exf  nation.  When  large  proteins 
and  other  aggregated  molecules  are  dissolved  in  water,  the 
volume  of  the  solution  is  usually  larger  than  the  volume 
of  solvent.  This  "volume  factor"  has  been  found  to  be  94 
percent  for  normal  plasma;  that  is,  in  Table  II,  the  28.50 
moles  of  HjO,  which  is  516  cc  at  38°C,  represent  516/ 

0.94  -  549  cc  of  plasma  per  liter  of  blood.  The  red  ceil 
volume  is,  thus,  451  cc  per  liter  of  blood  and  the  hemato¬ 
crit  is  45.1  percent.  Since,  in  the  red  cells,  there  arc 
only  18.0  moles  or  326  cc  of  H^O,  the  volume  factor  is 
72.3  percent,  which  Includes  the  volume  occupied  by  the 


stroma. 
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3.  DEVELOPMENT  OF  THE  BLOOD  MODEL 

GENERAL  CONSIDERATIONS 

The  data  of  the  previous  section  enable  the  Simula* 
tion  of  the  average  arterial  _>locd  of  Table  II,  which 
represents  a  liter  of  bio  .  the  steady,  resting  state. 

We  assume  that  metabolism  continues  at  a  steady  rate 
and,  because  metabolism  in  the  red  cells  is  principally 
anaerobic,  gives  a  net  production  of  H+  ions  and  lactate 
ion  along  with  other  minor  products.  We  also  assume  a 
complex  but  undefined  activity  (referred  to  as  the  "sodium 
pump")  for  the  red  cell  membrane.  One  result  is  that  all 
mobile  ions,  particularly  Ne‘  and  K  ,  acquire  across  this 
membrane  steady  concentration  gradients  different  from 
those  without  the  active  membrane.  Also,  the  system  has 
a  steady  "membrane  potent! al"--roughiy,  that  electrical 
potential  which  is  measured  by  suitable  electrodes  suit¬ 
ably  placed  (Chap.  Ill  analyses  this  definition). 
steady  pH  exists  in  the  plasma  and,  because  of  the  activity 
of  the  cell,  a  different,  lower  pH  holds  in  the  "interior" 
of  the  red  cell. 

The  structure  of  the  red  ceil  is  simpler  than  other 
body  cell*  since  the  red  ceil  is  not  nucleated  and  contains 
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but  one  principal  protein.  In  the  model,  the  interiors 
of  all  red  cells  in  the  liter  of  blood  are  taken  together 
as  a  single  !,compartment uniformly  mixed  and  separated 
from  the  plasma  by  an  ideal  semi-permeable,  though  active, 
membrane. 

Conceptually,  thin,  two  distinct  compartments,  with 
chemical  analyses  given  in  Table  II,  communicate  with 
each  other  through  the  idealized  membrane.  Proteins  and 
certain  other  large  molecules  (e.g.,  complex  phosphates) 
cannot  penetrate,  but  other  ions  and  molecules  pass 
through  regardless  of  the  mechanism  or  of  the  time  required. 
The  membrane  is  flexible  and  will  not  support  a  pressure 
gradient,  so  the  "red  cell"  may  change  volume  by  as  much 
as  5  percent  between  arterial  and  venous  blood.  Also, 
the  temperature  is  supposed  uniform  throughout. 

STATEMENT  OF  THE  PRELIMINARY  PROBLEM 

Given  these  basic  hypotheses  and  assumptions  (as  well 
as  further  assumptions  carefully  defined  as  i.quired),  the 
essential  problem  of  this  simulation  is  to  create  Table  II 
from  Table  I  by  means  of  a  mathematical  model--that  is. 
given  the  total  composition  of  a  liter  of  whole  blood  and 
using  the  general  mathematical  program,  distribute  these 
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"input"  substances  into  the  compartment’s  plasma  and  red 
cell  so  that  the  "output"  duplicates  Table  II. 

FIRST  APPROXIMATIONS 

The  construction  of  a  model  of  arterial  blood  could 
begin  with  the  simulation  of  ordinary  "soda  water."  a 
two-phase  system  of  gas  and  liquid  with  both  inter-  and 
intraphase  reactions. *  The  gases,  0^,  C^,  ^ ,  have 
determined  solubility  constants  at  a  given  temperature 
and  pressure,  and  the  solvent  water  has  known  vapor  pres¬ 
sure.  One  could,  therefore,  compute  the  concentrations 
of  each  in  the  liquid  phase  as  well  as  the  partial  pres¬ 
sures  of  IH^O  vapor  in  the  gaseous  phase.  In  the  liquid 
phase,  the  carbonate  system  reaction  constants  are  well- 
documented,*  so  with  fixed  partial  pressures  of  the  gases 

—  at 

one  could  compute  the  concentrations  of  HCO ^ ,  CO^,  H  , 

OK  ,  and  1^0  with  arbitrary  accuracy. 

If  the  gas  volume  is  finite,  the  formation  of  the 
species  in  solution  will,  by  conservation  of  mass,  alter 
the  gas  composition.  Therefore,  the  preliminary  model 
uses  a  very  large  gas  phase,  which  also  gives  better 

♦ 

Reference  ?  further  details  this  model. 

t 

See,  for  example  F.dsrll  and  Wyman,  Ref.  8,  Chap.  10. 
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control  of  independent  variables.  Using,  say,  1000  moles 
of  gas  phase  at  a  determined  composition  versus  one  liter 
of  liquid  phase  causes  the  gas  phase  to  remain  almost 
constant,  thus  holding  the  gases  in  solution  at  constant 
concentration.  The  gas  composition  can  be  varied,  of 
course,  as  required.^ 

MODEL  A- -PRIMITIVE  SYSTEM.  PLASMA 

The  first  model,  Model  A,  only  slightly  enlarges  the 
soda  water  model  which  consists  of  the  gases  0 CC^, 

H^O,  the  solvent  ^G,  and  two  phases,  gas  and  liquid. 
Adding  more  species  to  the  liquid  phase  without  adding 
intraphase  reactions  will  change  the  results  for  soda 
water  only  slightly.  Thus,  adding  Na+,  Cl  ,  K+,  HP0~, 

"H"  |  |  ss 

Ca  ,  Mg  ,  SO^,  urea,  glucose,  lactic  acid,  and  protein 
up  to  the  approximate  concentrations  of  these  substances 
in  plasma,  changes  the  concei  rations  of  dissolved  gases 
only  slightly,  although  the  pH  may  require  adjustment 
with  NaOH  or  HC1.  The  gas  concentrations  do  change  be¬ 
cause  the  solubility  constant  for  the  gases  differs  in 

+Frevious  studies--Ref3 .  1  through  5--describe  the 
computation  procedure  employed  here. 
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plasma  and  in  water,  and  because  the  actual  number  of 
moles  of  liquid  phase,  x,  has  increased.  No  new  intra¬ 
phase  reactions  are  introduced  yet.  Later  the  protein 
will  buffer  H+  ion  and  bind  Ca++,  Mg"1^,  and  Cl". 

Here,  the  protein  makes  the  solution  non- ideal  and 
increases  the  volume  of  plasma  over  the  water  volume 
by  6  percent.  Otherwise,  in  the  normal  viable  range, 
this  two-phase  preliminary  system  can  be  assumed  to  be 
well-behaved.  Again,  concentration  curves  could  be 
computed  for  this  simple  system,  but  it  is  not  yet  of 
much  interest. 

MODEL  B—OSMOLARITY 

A  second  model,  Model  B,  has  a  third  compartment 
called  "red  cells"  and  several  additional  conditions. 
For  this  model,  the  total  amounts  of  substances  added 
to  the  total  system  will  be  exactly  those  of  Table  I, 
arterial  blood.  The  problem  is  to  determine  how  much 
of  each  of  these  substances  is  distributed  into  the 
red  cell  compartment  under  various  conditions.  Such 
a  determination  is  a  solution  of  the  problem  under 


those  conditions. 
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First,  for  Model  b,  let  the  membrane  between  plasma 
and  red  cells  be  permeable  to  all  species--including 
protein.  In  this  case,  the  mathematical  solution  is  not 
unique,  for  the  membrane  could  divide  the  two  compartments 
arbitrarily,  and  the  mixed  composition  of  the  two  com¬ 
partments  would  be  identical  and  uniform  except  for  Xp^ 
and  the  total  number  of  moles  of  all  species  in 

each. 

Second,  Model  3  a)  supposes  the  membrane  impermeable 
to  proteins  only,  b)  ignores  the  electrical  charge  on  the 
ions  and  proteins,  and  c)  supposes  (as  do  subsequent 
models  in  Chap.  I)  that  the  gases  have  equal  solubility 
in  red  cells  and  plasma. 

With  the  protein  distributed  and  fixed  in  plasma  and 
red  cells,  and  with  the  membrane  flexible  so  that  hydro¬ 
static  pressure  gradients  cannot  arise,  all  species  dis¬ 
tribute  between  the  compartments  in  proportion  to  the 
amounts  of  impermeable  proteins  in  each,  a  result  of  the 
elementary  "osmotic"  phenomenon:  i.e.,  each  substance, 
including  water,  will  distribute  no  that  it  has  the  same 
concentration  on  each  side  of  the  membrane,  independent ly 
of  the  other  species.  In  addition,  since  the  proteins 
are  impermeable,  all  other  species  must  move;  i.e,,  ‘:he 


i 
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membrane  must  move,  so  that  the  proteins  on  each  side  are 
equally  dilute.  Thus,  the  total  moles  in  each  compartment 
have  the  ratio 

xpl/Xrc  *  S-706xl0_y(10. 568>  l(f 3  +  2.2725xl(f3)  -  0.0678 

where  the  numbers  on  the  right  are  the  moles  of  protein 
on  each  side  of  the  membrane.  Table  III,  a  reproduction 
of  the  computer  results  for  Model  B,  illustrates  these 
comments . 

MODEL  C-- ELECTRICAL  NEUTRALITY 

Model  C  Involves  an  hypothetical  experiment.  Suppose 
that,  in  addition  to  the  above  conditions,  each  protein 
now  has  an  electrical  charge  of  -1  per  molecule,  and  the 
other  ions  have  their  normal  chemical  valence.  Suppose, 
further,  that  the  sum  of  the  charged  input  substances  is 
such  that  the  whole  is  neutrally  charged,  Now  apply  th? 
mathematical  condition  that  each  compartment,  plasma  and 
red  cells,  be  neutrally  charged,  i.e.,  that  the  sum  over 
all  charged  species  in  each  compartment  is  zero.  Under 
these  conditions,  no  change  at  all  will  occur  in  the 
computed  distribution  of  the  species  and  the  results  will 
be  identical  to  Model  B.  That  is,  under  the  conditions 
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Table  III 

ELEMENTARY  OSMOLARITY- -MODEL  B 
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for  B  but  with  all  protein  uniformly  charged,  ail  sub¬ 
stances  are  already  distributed  on  both  sides  of  the 
membrane  such  that  each  compartment  is  electrically 
neutral. 

MuDEL  D--GIBBS-DONNAN  AND  HAMBURGER  SHIFT 

Model  D  supposes  that  the  average  charge  of  the 
protein  molecules  in  plasma  only  is  -10  per  molecule.  The 
moles  of  input  substances  are  identical  to  Table  I,  and 
the  average  charges  on  hemoglobin  and  the  miscellaneous 
fixed  species  in  red  cells  are  such  that  the  total  inputs 
are  neutrally  charged.  Appl>  ng  the  neutral  charge  re¬ 
straint  now  to  plasma  and  to  red  cels,  if  the  average 
charge  per  molecule  of  the  red  cell  protein  happens  also 
to  be  -10,  the  results  are  agaii  as  in  Model  B  or  C.  Any 
other  charge  on  the  average,  red  cell,  protein  molecule 
causes  ’Ifferent  results. 

The  average  charge  for  all  red  cell  proteins  in 
Table  I  is  -3.5,  and  a  comparison  of  this  with  -10  for 
the  average  charge  cf  plasma  protein  indicates  that  posi¬ 
tive  ions  should  be  displaced  into  plasma  and  the  permeable 
negative  ions  displaced  or  repelled  into  red  ceils.  Table 
IV,  a  reproduction  of  the  computer  results  for  Model  D, 


r 
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Table  IV 

GIBBS-DONNAN  PHENOMENON- -MODEL  D 
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shows  this  actually  occurs.  Table  IV  demonstrates  that 
Na+,  say,  is  now  predominantly  in  plasma,  and  Cl  is  pre¬ 
dominantly  in  red  cells.  The  total  mole  ratio  of  plasma 
and  red  ce’ 1  compartments,  hflS  c^ange^  t0  0.1455. 

Note  that  this  movement  of  the  positive  and  negative 
ions  is  anomalous.  The  directions  of  movement  of  the  Cl 
ion  and  K+  ion  here  are  opposite  to  the  standard  blood 
distribution.  Thus,  the  zero-charge  restraint  plus  the 
impermeable  fixed-protein  distribution  alone  are  in¬ 
adequate  to  establish  the  desired  (Table  II)  electrolyte 
distribution.  This  observation  is  significant  because 
Model  D  contains  both  of  the  essential  elements  for  the 
classical  Gi'bbs-Donnan  equilibrium  [9],  i.e.,  the  fixed- 
protein  and  the  charge  restraint. 

The  Gibbs-Donnan  phenomenon  is  usually  explained 
with  an  hypothetical  system  where  electrically  charged 
protein  is  on  only  one  side  of  a  rigid  membrane  in  a 
water  solution  containing  other  electrolytes.  The  zero- 
charge  constraint  is  then  applied  to  a  compartment  on 
one  side  of  the  membrane.  Substances  move  across  the 
rigid  membrane,  in  response  to  both  the  osmotic  effect  of 
the  protein  and  the  change  condition,  until  a  hydrostatic 
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pressure  builds  up  on  the  protein  side  sufficient  to  stop 
further  flow.  Simultaneously,  the  electrolytes  rearrange 
to  satisfy  the  charge  condition.  The  hydrostatic  pres* 
sure  thus  acquired  by  the  system  is  an  empirical  measure 
of  the  "osmotic  pressure"  of  the  system. 

In  the  present  case,  the  membrane  is  flexible  so 
hydrostatic  pressure  cannot  build  up--the  complementary 
compartments  simply  change  size.  Replacing  the  hydro¬ 
static  pressure  which  balances  the  flow,  a  second  protein 
is  dissolved  in  the  second  compartment.  Osmotic  re¬ 
distribution  will  cease  when  the  two  proteins  are  equally 
dilute,  as  in  Model  B.  The  charge  restraint  simultan¬ 
eously  redistributes  the  electrolyte,  as  in  Models  C  and  D. 
The  Gibbs-Donnan  phenomenon,  however,  is  basically  the 
same,  and  the  above  results  show  it  is  insufficient  in 
itself  to  establish  the  desired  (Table  II)  distribution 
of  electrolytes  for  arterial  blood. 

The  Gibbs-Donnan  equilibrium  phenomenon  does,  how¬ 
ever,  contribute  an  important  characteristic  to  the  blood 
model--namely ,  the  so-called  "chloride"  or  "Hamourger 
shift"  which  takes  place  reversibly  at  the  lung  surface 
and  in  the  body  tissue  capillary  beds  [l  .  At  the  body 
tissue  level,  CO^  dissolves  in  plasma  and  diffuses  into 
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the  red  cell,  where  the  en.  yme  carbonic  anhydrase  com- 
bines  it  with  a  molecule  of  water  to  become  l^CO^.  The 
consequent  ions  dCO^  and  H+  are,  of  courte,  subject  to 
the  charge  constraint  of  the  Gibbs-Donnan  phenomenon, 
just  as  any  other  ion.  Conventionally,  the  charge  con¬ 
straint  affects  the  concentration  of:  HCO^  exactly  as  it 
does  Cl  ,  yielding 


CHCOl]  .  tCl']  . 

- i-£±  .  k - —  (1) 

rHc°;  jRC  [cnRC 


where  square  brackets  indicate  concentration.  In  this 
case  k  ■  1;  Model  G  below  introduces  a  difference  in 
solubility  for  CO2  between  plasma  and  red  cells,  so  k  ft  1. 
As  the  concentration  of  bicarbonate  increases  in  the  red 
cell,  the  HCO^  moves  into  the  plasma  under  a  diffusion 
gradient.  This  movement  o?  HCO^  would  leave  a  net  posi¬ 
tive  charge  in  the  red  cell  except  that  Cl  ,  as  well  as 
other  anions,  then  shifts  into  the  red  cell  milieu  (the 
chloride  shift).  The  solvent  species  must  also  shift 
to  maintain  Eq .  (1).  Exactly  the  opposite  set  of  re¬ 
actions  occurs  at  the  lung  surface. 
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Two  important  tentative  assumptions  derive  from  this 
discussion.  First,  since  Eq.  (1)  is  accepted  in  the 
literature  [101  and  verified  empirically  [6,7 neither 
Cl  nor  HCO^  is  likely  to  be  subject  to  active  transport 
(see  below,  Model  E)  at  or  in  the  cell  membrane.  If  the 
concentration  ratios  of  either  were  effected,  Eq.  (1) 
could  not  hold  over  the  range  of  the  Hamburger  shift  un¬ 
less  both  were  affected  exactly  equally  by  an  active 
transport  meenanism,  which  is  unlikely. 

The  second  assumption  is  that  Eq.  (1)  may  be  made 
much  stronger.  Indeed,  any  ion  not  subject  to  active 
transport  must  have  the  same  chemical  potential  gradient 
given  by  the  zero-charge  constraint,  i.e., 


[Cl’l 


PI 


k 


rHC0']pl 


r  k+ /  r so”  i 

h  RC  i  -  u4  PI 


h 


[Ci‘ 


RC 


[HCO 


3  RC 


[H 


(2) 


PI 


tso.i 


4  RC 


where  the  hydrogen  ratio  is  inverted  because  of  its  sign 
and  the  SO^  ratio  is  reduced  because  of  its  valence. f 
The  dots  indicate  similar  ratios  for  other  ions.  For 

4-  + 

example,  in  the  present  Model  P.  Na  and  K  are  not  subject 
*See  Ref.  8,  p.  227. 
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to  active  transport  (see  Model  E)  and  therefore  they,  too, 
satisfy  Eq,  (2),  as  shown  in  Table  IV. 

MODEL  E— THE  SODIUM  PUMP 

Model  E,  again  a  hypothetical  experiment,  ignort 
the  electrical  charge  of  the  species  in  order  to  simplify 
a  first  simulation  of  the  sodium  pump.  Because  the  species 
subject  to  the  sodium  pump  have  electrical  charge,  the 
effects  of  the  pump  and  the  effects  of  the  ?,ero-charge 
restraint  applied  in  D  would  otherwise  interact.  In 
order,  then,  to  isolate  the  function  of  the  sodium  pump, 
the  mathematical  serc-charge  restraint  of  the  Gibbs- 
Donnan  phenomena  is  te-maorariiy  suspended.  All  other 
condi. tions--semipermeable ,  flexible  membrane,  etc. --still 
apply. 

"Sodium  pump"  is  a  generic  name  for  the  phenomena 
whereby  living  cells,  via  a  biochemical  and/or  biophysical 
mechanism  within  or  at  the  surface  of  the  membrane  of 
those  cells,  selectively  excrete  or  absorb  cations.  The 
sodium  pump  is  a  subclass  of  the  more  general  "active 
transport"  processes,  the  exact  mcchanir.ms  for  which  are 
still  generally  obscure  "11,12\ 
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The  present  Model  E  concerns  not  the  biochemical  or 
biophysical  mechanism  whereby  the  sodium  pump  is  accom¬ 
plished,  but  only  the  net  result.  As  may  be  seen  in 
Table  II,  the  sodium  pump  acting  in  the  erythrocyte  mem¬ 


branes  causes  most  of  the  Na  and  Ca  ions  to  concentrate 


in  plasma,  and  most  of  the  K  and  Mg  ions  to  concentrate, 
in  the  interiors  of  the  red  cells.  The  above  examples 
show  clearly  that  the  Gibbs-Donnan  forcing  functions, 
osmotic  effects,  and  zero-charge  condition  cannot  account 
for  the  disposition  of  cations  in  Table  II.  Additional 
work  on  these  cations  is  necessary. 

In  order  to  quantize  the  sodium  pump  phenomena,  we 
assume  that  the  species  ^apu^SI^  and  the  species  Na^ED  CELL 
are  not  identical.  It  is  convenient  to  say,  instead,  that 


Napv :  n4  • 


(3) 


is  a  chemical  reaction  converting  one  species  to  the 
other,  an  interphase  reaction,  even  though  Eq.  (3)  is 
not  a  complete  reaction  in  the  thermodynamic  sense, 
further,  the  increment  of  free  energy,  AF,  is  proportional 
to  the  work  done  on  the  substance  during  conversion,  al¬ 
though  the  actual  work  in  vivo  is  difficult  to  measure 
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and  may  differ  considerably  from  this  theoretical  minimum 
v^lue  required  to  maintain  a  gradient.  Thus,  though  we 
do  not  represent  the  membrane  cr  the  mechanism,  we  postu¬ 
late  a  free  energy  parameter  proportional  to  the  dif¬ 
ference  in  concentrations  or  mole  fractions  of  the  two 
species  in  Eq.  (3)  at  steady  state. 

Tables  V  and  VI  respectively  list  the  input  data 
and  reproduce  the  computer  results  for  Model  E.  Table  V 
lists  the  species  to  be  expected  in  the  three  compartments, 
their  respective  free  energy  parameters,  and  the  input 
components  from  which  the  output  species  originate,  as  in 
a  chemical  reaction. 

The  "free  energy  parameters"  in  Table  V  have  various 
appropriate  meanings.  The  species  has  a  free  energy 
parameter  -10. 94  in  the  "Air  Out"  compartment,  and  zero 
in  "Plasma."  This  number  is  the  log  bsse  e  of  the  solu¬ 
bility  constant  in  mole  fractions  for  09  in  plasma, 

0.02089  co/cc/af.m.  The  species  HCO^  in  plasma  has  the 
free  energy  parameter  -21.35,  which  is  the  log  base  e 

for  the  mass  action  constant  in  mole  fractions  for  the 
♦ 

t  faction 

OH’  +  CO.,  -  HCC~  ,  AF  (4) 


7 - 

See  Ref.  2 ,  p.  15 . 
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Table  V 

INPUT  DATA- -ELEMENTARY  MODEL 


Free 


Expected 

Energy 

Input 

Species 

Parameter 

Components 

*1*  OUT 

0? 

-10.«40000 

1.000  02 

CO? 

-t. 690000 

1.000  C02 

H2 

-11. 620000 

1.000  N2 

M20 

- 16 . 600000 

1 .000  H* 

1.000  OH- 

Pl*S*« 

02 

0. 

1.000  02 

cn? 

0, 

1.000  C02 

N? 

'J. 

1.000  N? 

H» 

0. 

1.000  «♦ 

PH- 

0. 

1.000  OM- 

Cl- 

c. 

1.000  Cl- 

Hk* 

0. 

1.000  **♦ 

K  ♦ 

0. 

1.000  K* 

C  **♦ 

0. 

1.000  C*44 

MG** 

0. 

1.000  HG*» 

S04* 

0. 

1.000  S04- 

HP04  ■ 

n. 

l.OGJ  MPQ4* 

U*M 

0. 

i.ooo  mu* 

r.iuco* 

0. 

1.000  GIUCOS 

l *CTtC 

0. 

1.000  L*CT!C 

NH4* 

0. 

1.000  NM4* 

►X01- 

-21 . 160000 

:.ooc  cu2 

1.000  OM- 

h?CP1 

-32.840000 

1.000  C02 

l  .000  H* 

CO'* 

6.260000 

1,000  CO? 

-1.000  H* 

W?0 

-1«J.  180000 

1.000  H* 

1.000  OH- 

PISCPP 

0. 

1.000  *!VC*L 

■  to  cuts 

p? 

-0. 

1.00<-  02 

CO? 

0. 

1.000  CO? 

N? 

-c. 

1.000  H! 

M  ♦ 

0. 

1.000  H* 

OM- 

0. 

1.000  ()N- 

Ct- 

0. 

i.oeo  u- 

Hk* 

2.  1  'Mir, 

1.000  Hk » 

■  ♦ 

-  .  >  %  1  ! 

1 .000  *♦ 

f 

i  .  .’4  1  It  o 

1.000  <.*♦♦ 

•  G  »  • 

-‘i.kSt  in  i 

1 .000  ac«» 

404* 

. 

l.OOC  104. 

**P04« 

-P. 

I.ooo  M*04- 

* 

0. 

1 .000  CrtM* 

GtOCOS 

n. 

1.000  GilXOI 

licnc 

0. 

i.ooo  i*cm 

vw  4  • 

0. 

1.000  NM4* 

HC01- 

-  ?S  .  160000 

1.000  CO? 

1.000  OH- 

**2CP  1 

- >2 . *40000 

1.000  CO? 

l  .  -  "00  H* 

CO  l  * 

6 . 240000 

1.000  CO? 

-  1  .  OOO  **1 

**?fl 

18.  180000 

1 . 000  H* 

1  .000  OH- 

c . 

I.ooo  «ISCM 

HPi 

0. 

i.OOO  H*4 

1.000  OH- 
I .000  OW- 


1.000  U*»- 

l .  OOC  P**- 


JL 


-29- 


Table  VI 

SODIUM  PUMP- -MODEL  E 


«Kt"IN**Y  *1000  TOOK  TO  IMMG7  t 00*110*  HIMO 


•IK  OUf 

ataiM 

K(0  Citll 

Mia 

*.**11*2  0/ 

2.**)l?f  01 

1. *•?*)€  01 

OH 

7.111*77  00 

Minn  oo 

0? 

TOifs  i.mort  0/ 

*,2«t)ti~oi 

i.Mom-oi 

"M»C  1.11*707-0! 

7. 111227-0* 

2.111221-0* 

CO/ 

•*01  f  S  1.H0//I  01 

4.1*0*27-01 

1.  7*4*77  -0* 

"M*C  1.2**777-0/ 

2.*o7***~os 

2. 107**7-0* 

HI 

"01 7  i  7.110007  0/ 

7.01*7)7-01 

1.17**17-01 

"f*»r  7.11*177-01 

7.1*1)27-0* 

7.1*1)27-0* 

H/0 

"Ot  f  1  6.10/717  01 

.’.Mom  oi 

1. *7**07  01 

"f  *»C  6.107**7-02 

*.*11127-01 

*.**11/7-01 

Ht 

"0171  -0. 

M*tlli-0t 

1.* 1*007-0* 

"7**C  -0. 

•.111077-10 

a.  111077-10 

OH- 

mil  -0. 

MIHIMI 

1. **7*17-07 

■Y**C  -0. 

*.1172)7-0* 

*.1172)7-0* 

Cl- 

"OUS  -0. 

i.*o7oa*-oz 

1.1*7*27-02 

"7*»C  -0. 

1.70*2*7-0) 

1.70*2*7-0) 

"*♦ 

"OUS  -0. 

7. *172*7-0? 

*.**7217-01 

"7*»C  -0. 

7.*07**7-0) 

1.7*1117-0* 

"OUS  -0. 

1.00*701-0) 

1.20*1 17-02 

"FK*C  -0. 

1. 11*2)1-01 

2.111*07-01 

C»»» 

mm  -o. 

i.*m.'t-oi 

1.1)17*7-0* 

"7**C  -0. 

l.ai  1707-0* 

«c«* 

»ous  -o. 

7.iii*ii-oi 

a  .***i*7  -o* 

"YK*C  -0. 

7. 7*/l*|-01 

1.  >*)«*•  -01 

10»« 

•OUS  -0. 

l.O/llli-Ol 

7.77*1*7-0* 

•»#••<■.  -0. 

1. 121001-0* 

1.171007  -01 

HOO*' 

"01  IS  -0. 

a. 2 looi i-o* 

*  .0*  ***7  -0* 

"»**C  -0. 

1.0111*1-0* 

1.0)11*7-0) 

UMI 

TOUS  -0. 

l.*0«)0»-0l 

l. 1)1707 -0) 

«»**c  -a. 

o.  7o**u-oi 

». 70**17 -0) 

UUCOl 

•oi f  s  -o. 

/. io*a*i-o i 

1 .111171-0) 

•*1*C  -0, 

7. *277*1-01 

F.K777M-01 

l»C  TIC 

•f'l  » %  -  o . 

1. 1 70011-01 

*.•7*177  -01 

■V**l  -0. 

1. 110*11-0) 

*. I)»*«7 -01 

•Ol  11  -0. 

1. 0O» It* -01 

).*•*•*« -0* 

mu  HI  -n. 

1 • *11)77 -Cl 

t . *1)177 -Cl 

atm- 

•WM  -P. 

t . t !**•! -  07 

•  .111IH  -01 

•MU  -0. 

1 .  1071*7-0) 

*.  K  7**7  -0* 

"7C0> 

»>*M  -0. 

• .  7  101 7f  -C  7 

*.*01117  -01 

•  *C  -  P . 

1 . 1711*1 -0* 

1  ■  *7*1*7  -0* 

t  01» 

mi\ 

1 .111717-01 

1.0171*7-01 

will  -0. 

1.  11)707  -  07 

i.  )*i;o*-oi 

•  IU»' 

•01  *  1  -  '•> . 

7.  70*007  -01 

1 .01**0*  *7 

<MK*C  -B. 

1.77***7-01 

1.  H  *027  -0* 

»•* 

•oi  a  i 

-8. 

7.7  77*07  -0 1 

•tiu 

-P. 

1.1*11*1-0* 
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plus 


giving 


H20  +  C02  -  H2C03  -  H+  +  HCO* 


OH*  +  H+  ~  H20 


(HCO‘) 

(oh*)(co2) 


=  K  . 


Using  the  free  energy  parameters  for  Na£  and  Na+ 

r '  KC 

in  Table  V,  we  may  compute  from  Table  VI 


n  (N«;1y(Na; 


,)  -  In  (2.908*10  73. 243*10‘4)  «  2.1934 


That  is,  the  cations  have  acquired  the  gradients  given  by 
their  respective  free  energy  parameters,  c^.  proportional 
to  AF  of  Eq  ■  (3).  Taking  their  valence  into  cons ioerat ion , 
the  ^ther  cations  satisfy  similar  eqiations.  References  1 
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and  3  develop  the  detailed  mathematical  background  for 
the  procedures  and  computations  involved  here. 

Note  that  the  distribution  of  species  in  Table  VI 
is  not  yet  that  of  Table  II.  From  Table  VI,  because  of 
the  new  disposition  of  ♦ he  solute  resulting  from  the  cation 
pumps,  the  ratio  of  Xp^/Xpc  is  now  26.95/19.88  -  1.355. 

The  pH  *  7.34  is  the  same  for  each  compartment  since  H+ 
ion,  assumed  to  be  not  subject  to  an  active  pump,  is 
distributed  uniformly  throughout.  Model  F  again  applies 
the  charge  constraint,  causing  the  H'  ion  to  shift  nearer 
to  that  of  Table  II. 

MODEL  F-- SIMULTANEOUS  CONDITIONS 

Model  F  finally  applies  the  conditions  of  the  previous 
examples  simultaneously,  i.e.,  the  cation  pumps,  the  zero- 
charge  restraint,  and  the  impermeable  protein  with  average 
charge  -10  per  molecule  on  the  plasma  proteins. 

Ine  input  data  for  this  model  are,  therefore,  es¬ 
sentially  the  same  as  Table  V,  except  for  the  equation 
which  computes  and  imposes  the  charge  constraint.  Table 
VIX  reproduces  the  computer  results  for  Model  F. 

An  anomalous  shift  of  electrolytes  has  Important  re¬ 
sults.  Compared  with  Table  III  the  charge  equation 
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Table  VII 

NONLINEARITY— MODEL  F 


*Atll«!7*A*T  ftlOOO  ncrtl  A07M  CMAAGf  AND  N*  At m* 


» t •  out 

At  ASHA 

A*.  lilt 

A.AAtIA*  04 

l.tIWII  01 

i.7»rot*  01 

*M 

-0. 

7 ,  4000 If  00 

7 .4)4*4*  CO 

0? 

*oif>  i.mo/i  oi 
•fa  ac  t.itom-ot 

». ioom-0* 

4.  ) >444? -04 

4.  l4t*4t-04 
4.  ))4?4«-04 

eo; 

"CHS  4.44 04/|  ot 
•»*AC  4.44477*-04 

7  •  0404AF-04 
4.40  7«*0*-0*. 

4 .44  7 1 11-04 
4 • 40794* -0) 

H2 

•01*  S  7.440UO*  0? 

•rite  r.AAAA/f'Oi 

7. 1  AAF.lt -04 
7.44))4*-l>6 

t.  74474* -04 
7.447)4*  -0* 

m*o 

•OlfS  6.10/Mt  01 
•  **AC  7>.  107947-0  4 

4. *9441*  Ot 
4. 94444* -Ot 

1.74 AOAF  Ot 
A. 44444* -0 t 

M* 

•OlfS  -0, 

•F4  4C  -0. 

4.0911 |  F-CA 
7.|744lf-tO 

1 .14444* -00 
l  .04*9)* -04 

OH- 

•nits  -o. 

"FAAC  -o. 

1. 160401-07 

1 .OAAOIf-OA 

1. 710407-07 
7.414*1* -04 

Cl  - 

•Ot  F  A  -<». 

•FAAC  -0. 

4.44944*-04 

l.F40«IF-01 

4 .  )4A  )  A*  -0? 

1  .  )4  7,4 1*  -0) 

Nf  ♦ 

•01  Ft  -0. 

•FA  AC  -0. 

7. 71A14F-04 
4. 447401-0) 

r.t)«)li-0) 

4.71  A44*  -04 

*• 

•Ot  ft  -0  . 

•F  A  AC  -L). 

MtOMt-O) 

A.4AOMF-04 

4 . 4)  14  )t  -04 
4.40404* -0) 

(Am 

•nits  ->}. 

•FAAC  -0. 

t  .4|4»0*-0I 
4.44444* -04 

l.  44  70)*  -04 
t ."4000* -04 

nc  ♦  * 

•••It  F  A  -0. 

•FAAC  -0. 

4.4»'44*-04 

I  .At  )4)f -04 

i.0tl74*-0) 
4.  1  )0t  tl -04 

•Cl *  A  -0. 

•  FAAC  -0. 

4.040441  -04 
» .40447* -04 

t.  14474*  '4 

4. 4*004* -0* 

■wit  1  A  -0. 

•FAAC  -". 

t  .  t  1 4101 -0  t 
1.1414)4  -04 

t,  t4»»4|  -0* 

1  .  7*4*1*  -04 

* 

«ot»4  -0. 

•  FA  AC  -0. 

i.mttt-oi 

A. 70»A|* -04 

t  .1*40  7*  -0) 
4.  70*41*  -»4 

04  IXU\ 

"Oil  7  -.1. 

•FAAC 

4.4#4l4»-t« 
7.  A.1  71*1  -04 

t  .  **)••«  -O 

7  .  *4  7  7M  -04 

IUYK 

•am  -0. 

•FAAC  -0. 

1  .4  1  ItOF  St 

4  .  *44  1;  |  -134 

4.44404*  -04 
I.  t  7  044*  -04 

H**%  • 

•Alt 
•»••'  -n. 

4  .  A .144*1  -04 

1  .  4  747»7  -04 

4  .  ‘70  4  1 4*  04 

4,  II  >4  7*  -04 

«C  D%- 

•Hit  -o. 

•JAAC  -5. 

1 . 444441 -04 
* . *44441 -04 

4.404  14*  -0* 
t.  *4341*  -04 

••i*  ft 

Mill  -rt. 

A . 444041  jt 
t  .  4444*1  -  04 

4, "l'4*4t  -fl 
1.  4/4444  -04 

tr  »* 

■rt i %  -s. 

•*  A  At  - 

4 . 04  7*4* -34 
4.  Ft  4*4*  -  0  7 

4.  74  >«|*  -0 A 
t.  ft)  If*  -07 

•  !V(H 

•tit  ■?. 

•••«<  - 

A.  -S4 

; .  ***;et -*t 

1.344*0*  -•/ 
4. «**••«  -04 

-4  4 

•ot  it 

"•It 

•0. 

4.47/40*  41 

1  .  .*444  « 

OflA'tO 
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effect  of  Table  IV  drives  a  small  amount  of  Na  Into  the 
plasma  compartment.  On  the  other  hand,  compared  with 
Table  VI,  the  a ame  charge  equation  effect  of  Table  VII 
drives  a  small  amount  of  Na+  into  the  red  cell  compart¬ 
ment.  That  is,  the  charge  equation  applied  to  Model  B 
results  in  a  different  direction  of  ion  movement  than 
would  the  same  charge  equation  applied  to  Model  E. 

Thus,  it  would  be  convenient  if  a  stress  or  forcing 
constraint  applied  to  the  system  always  gave  the  same  re¬ 
sults.  But  this  wish  is  naive.  The  result  of  a  stress 
is  clearly  a  function  of  the  system  to  which  it  is  applied. 

The  hematocrit  of  Table  VII  is  now  44  percent;  the 
pH  of  plasma  is  nearly  correct  at  7.4,  but  the  pH  of  red 
cells  is  still  alkaline  (compared  to  Table  II).  The  model 
requires  further  modification. 

MODEL  G-- HILL'S  EQUATION.  MYOGLOBIN 

Model  G,  the  final  model  in  this  preliminary  series, 
applies  three  small  modifications  tc  Model  F.  The  crude 
outlines  of  a  simulation  of  Table  II  are  already  apparent 
in  Model  F.  The  general  descriptive  character  is t ic«  of 
live  arterial  blood  .it  steady  state  have  been  deduced 
according  to  the  classical  principles  of  the  Moehemical 
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9tructure--i. e. ,  using  the  fixed  proteins  and  the  charge 
constraint  (to  obtain  the  Gibbs-Donnan  phenomena),  and 
the  sodium  pump  hypothesis.  At  this  point,  the  model 
resembles  that  of  Ref.  1  which  was  derived  from  a  dif¬ 
ferent  point  of  view. 

The  first  kind  of  addition  to  be  illustrated  con¬ 
cerns  additional  interphase  reactions--or  transport.  As 
an  example,  consider  the  HPO^  and  SO^  anions.  The  total 
of  these  ions,  as  acid- soluble  phosphate  and  sulfate, 
shown  in  Table  I  were  computed  from  Edelman  r 1 3 ^ ;  the 
distribution  shown  in  Table  II  was  estimated  from  Spector 
[7]  and  fv-elmar.  [13'.  The  anions  considered  here  do  not 
include  the  phosphates  and  sulfur  bound,  for  exampl  ,  in 
the  structure  of  the  nuc leoproteins  and  lipids 

In  Table  VII,  these  acid-soluble  anions,  in  response 
to  the  charge  equation,  have  a  greater  mole  fraction  in 
plasms  than  in  red  cells,  lust  as  the  chloride  ion  does, 
in  fact,  however,  most  of  f  r  e  H'  ,  in  whMe  blood  is  in 
the  red  cell  interior  rl\  and  we  presume  the  same  tor 
SO^  13’,  although  the  evidence  for  this  is  meager  and 
somewhat  contradictory . 

Ion  selectivity  ot  the  membrane  may  cause  these 
anU~s  *  o  appear  mainly  In  the  cell  interior;  HPC*  in 
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particular  may  become  involved  in  large  complexes  with 
low  permeability,  giving  *2  choice  of  methods  to  insure 
the  approximate  distribution:  An  intraphase,  red  ceil 
reaction  to  combine  the  HPO^  and  SO^  with  an  impermeable 
species;  or  interphase,  transp, rt- forcing  functions  similar 
to  the  scclum  puir.p.  Either  process  would  hold  these  red 
cell  aniens  in  a  predetermined  proport ion--a  proportion 
which  would  vary,  however,  with  the  number  of  other  anions 
present  because  of  the  charge  equation  and  Gibbs-Donnan 
effects . 

We  choose  for  this  example  the  interphase  pump  for 
distributing  the  HFO*  and  S0^\  Thus,  in  Table  VIII,  the 
free  energy  parumeteis ,  c^,  for  these  anions  differ  by 
2.0  between  plasma  and  red  cell.  As  in  Eq .  (6),  this 
causes  a  mole  fraction  gradient  between  the  two  compart¬ 
ments:  the  tnole  fraction  ratio  is  approximately  three, 

the  mole  fraction  being  greater  in  red  cells  (Table  IX). 
More  detail  will  be  added  in  later  models  where  the  phos¬ 
phate  also  ignites  in  a  buffering  reaction 

The  second  modification  gives  the  atmospheric  gases 
a  different  solubility  coefficient  in  the  red  cell  milieu 
than  in  plasma.  Using  the  data  from  Roughton  '14'.  Table 
VIII  shews  the  computed  differences  in  solubility  as 
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Table  VIII 

INPUT  DATA- -MODEL  G 


Expected 

Species 

Free 

Energy 

Parameter 

AIK  OUT 

02 

-10.940000 

C02 

-7  690000 

M2 

-11.520000 

M.O 

-36.600000 

PLASMA 

02 

0. 

CO  2 

0. 

C. 

H* 

0. 

OM- 

0. 

CL- 

0. 

0. 

K* 

0. 

CA** 

.1. 

HC'* 

0. 

S04« 

0. 

HP04* 

0. 

UREA 

0. 

GLUCOS 

0. 

LACTIC 

0. 

MH4  + 

\J  4 

HC01- 

-21.350000 

H2C03 

-  32 . 84U000 

CO  3* 

6. Z 60000 

H20 

-39.390000 

MISCPK 

0. 

RED  CELLS 

02 

-0.490000 

C02 

0. 

N? 

-0.500000 

H* 

0. 

OH— 

0. 

C'  - 

0. 

NA* 

2 . 193799 

K* 

-2.941575 

CA** 

Z.2517'70 

MG** 

-0.457703 

S04* 

-2.000000 

HP04» 

-2.000000 

UREA 

0. 

CLUCOS 

0. 

LACTIC 

O  ■ 

NH4* 

0. 

MC03- 

-21.490000 

H2C03 

-32.840000 

CO  3 

6.120000 

H20 

-39.390000 

MISCPR 

0. 

H84 

0. 

H6*02 

-  15.230000 

Input 

Species 


1 .000 

02 

1.000 

C02 

1.000 

M2 

1.000 

H* 

1.000  OH- 

1  .000 

□2 

!l .  000 

Ci>2 

i.000 

M2 

1.000 

H* 

-0. 

-0. 

1.000 

OH- 

-0. 

-0. 

1.000 

Cl- 

-0. 

-0. 

l.OOu 

MA* 

-0. 

-0. 

1.000 

*.♦ 

-0. 

-0. 

i.^OO 

CA** 

-0. 

-0. 

1.000 

MG** 

-0. 

-0. 

1.000 

S04» 

-0. 

•0. 

1.000 

MP04* 

-0. 

-0. 

’  .000 

UREA 

1.000 

GLUCOS 

1.000 

LACTIC 

-o. 

-c. 

1.000 

MH4* 

•  r  , 

-0. 

1.000 

C02 

1. 

000 

OH- 

-0. 

1.000 

C02 

1. 

000 

1.000  OM- 

1.000 

CO? 

_  P 
*  • 

000 

H ♦ 

1.000  OM- 

1.000 

H* 

l. 

OOil 

OH- 

1.000 

MISCPL 

-0. 

-0. 

1.000 

02 

1.000 

C02 

1.000 

N2 

1.000 

H* 

1.000 

OM- 

1.000 

:i- 

1.000 

N4* 

1.000 

K* 

1.000 

CA** 

l.OCO 

MG** 

1.000 

S04' 

1.0Cr 

HP04« 

1.000 

UREA 

1.000 

ucos 

1.000 

LACTIC 

l.nno 

MH4. 

1.000 

CO? 

1 

.000 

OH- 

1.000 

C02 

l 

.000 

H* 

1. 000 

OH- 

1.000 

CC2 

-1 

.000 

H* 

1.000 

OH- 

1.000 

H* 

l 

.000 

OH- 

1.000 

MISCRE 

l.Oi  0 

HB4 

1.000 

02 

1 

o 

6 

MB4 

1.000  »PLASM 
-1.000  •PLASM 
-1.000  (PLASM 
1.000  »PL ASM 
1.000  »PLASM 
2.000  »PL  SM 
2.000  (PLASM 
-2.000  »PL ASM 
-2.000  • PL ASM 


-1.000  (PLASM 
1.000  *PLASM 
-  1 .000  *P*. ASM 

-2.000  (PLASM 

-10.000  »PLASM 


A 


Table  XX 


COMPLETED  DISTRIBUTION- -MODEL  G 


preliminary  nmnn  model  sto  state  no  protete)  reactions  oei  and 


Aia  cut 

PLASPA 

RED  CELLS 

A-8AP 

9.99127E  02 

2 .  88  ?44E  01 

1.79S96E  01 

PM 

“0. 

7.  39264E  00 

7.19424E  00 

0? 

WILES 

t . 31 498E  0? 

6. 73934E-03 

6.842SSE-0S 

MFPAC 

1  .  1161  EE-OI 

?.  33409E-06 

I.8Q99  7F-0A 

CO? 

POLES 

S.26022E  01 

6.93300E-04 

4. 3/4  TOE -04 

pfrac 

S.26487F-02 

2 . 4080 1  E-OS 

2 • 40801E -OS 

N? 

ES 

I.S4000 E  02 

2.16367E-04 

2.2I882E-04 

P..C 

T.SA6S9E-01 

7.49S39E-  0 

1 . 23S4SE -OS 

M20 

POLES 

6.10?6>»E  01 

2.S7139E  Of 

l.78S9Sfc  O' 

PEAAC 

6.I0797F-0? 

9.94440E-01 

9.94440E-  . 

H* 

POLES 

-0. 

2.  10661E-08 

J.0690SEC8 

ME  0  AC 

-0. 

7.29S76E-I0 

l . 1 3206E -09 

OM- 

POLES 

-0. 

3.07724E-0T 

! .2121  IE -07 

PEPAC 

-0. 

1.06373E-08 

6.74909E-09 

ci- 

POLES 

•0. 

S. 7S290E-02 

2 . 262 lOE-O? 

PEPAC 

-0. 

1.98892E-03 

1.2595SE-03 

NA» 

POLES 

-n. 

7.644SSE-0? 

8.  37448E-03 

PFRAC 

-0. 

2.647S2E-03 

4.6629SE-04 

R» 

POLES 

-0. 

2.297SSE-03 

4.27S25E -02 

PEP  AC 

-0. 

7 . 9S704E-03 

2.38048E-03 

CA*» 

POL  f  S 

-0. 

1.37983E-03 

2 .2S IS 4E-04 

PEPAC 

-0. 

4.7787BE-OS 

1.25367E-0S 

PG  >♦ 

POLES 

-0. 

A  .  6  72  37E-04 

1 . 14S26E-03 

PF«  t 

-0. 

1 .61S17E-0S 

6.3  7bB8E -OS 

S04» 

"OL  f  S 

-  0. 

l • 84633F-04 

3.40347E-04 

PERAC 

-0. 

6*  39302E-Q6 

t  .89S0TE-0S 

HP04  * 

POLES 

-0. 

3.029S9E-04 

9.27041E-04 

P‘RAC. 

-0. 

l, 74188F-0S 

S.16I81E-0S 

UREA 

POLES 

-0. 

L.93713E-03 

1 .20487E-03 

PFRAC 

6. T0879E-0S 

6. 70879E-0S 

GLUCOS 

POL  E  S 

-0. 

2.26019E-03 

l  .4OS8IE-03 

t^RAC 

-0, 

7.R2764E-0S 

7.82  764E -05 

LACTIC 

POLES 

-0. 

1 .40R60E-03 

S.  74497E-04 

PEPAC 

-0. 

S.OSl i9E-0S 

3.  19883E  -OS 

NM4* 

POL  t  S 

-r . 

4.3B409F-04 

4. 30S92E-04 

PEPAC 

-0. 

I.S183SE-US 

7  .  397  SSE-OS 

HC03- 

POLE  S 

-  0. 

l . 3R678F-G? 

6.2J"3E-03 

PEPAC 

-0. 

4.80280F-04 

3  •  49RS9F  -04 

h?co  i 

POLES 

-0. 

9.88831 F-07 

6.  LS041E-0T 

PEPAC 

-0. 

3.424S9E-08 

3  .  4,'A  ->vt  -OB 

C03* 

POLES 

-0. 

1.9SU8F-03 

S. '  „9$8E -06 

PEPAC 

-n. 

6.  72284E-07 

».  .Jl  331  -07 

piscpp 

POLE  S 

-0. 

8. 706U0E-04 

3.7S0S0E 

PEPAC 

-0. 

3.01S12I -OS 

i  .  088  30E  -  °4 

MBS 

POL  E  5 

-0. 

•  0. 

3.  3S368E-04 

PFAAC 

-0. 

-0. 

1.867  34E -OS 

H6402 

MOL  f  S 

-0. 

-0. 

8. 734631: -03 

PFRAC 

-0. 

-0. 

4. R7442E -04 
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mod  if  led  ,  free  energy  parameters.  The  result  is  the 
solution  of  slightly  more  0^  and  ^  in  the  red  cell  milieu 
than  in  plasma.  Reference  2  computes  these  data. 

Finally,  this  preliminary  model  incorporates  an 
illustrative  reaction  of  hemoglobin  with  oxygen.  Hemo¬ 
globin  oxygenation  has  an  extensive  literature  (summarized 
in  Ref.  14)  and  is,  of  course,  still  the  subject  of  much 
research.  This  preliminary,  illustrative  model  incorporates 
one  of  the  earlier  theoretical  structures  of  the  oxygena¬ 
tion  reactions.  Studying  a  simpler  theory  and  its  short¬ 
comings  clarifies  some  of  the  details  of,  and  much  of  the 
motivation  for,  the  following  two  chapters. 

In  1910,  A.  V.  Hill  [15]  proposed  that  reaction  of  a 
solution  of  hemoglobin  with  oxygen  could  be  written 


Hb  +  nO  -  Hb(0„)  ,  K 

t  z  n 


(7) 


where  both  K  and  n  were  to  be  determined  from  empirical 
data  from  the  laboratory.  Using  the  Hill  equation,  as  in 
Ref.  14,  the  percent  saturation  of  hemoglobin  with  oxygen, 
y,  can  be  computed  from 


-X-  „  v* 
100  y 


V  ^ 

—M _ 

1  +  Kpn 


(8) 


where  p  is  the  oxygen  pressure  in  millimeters.  The 
saturation  can  be  measured  in  the  laboratory  at  various 
pC^  by  a  variety  of  means.  Comparing  the  laboratory 
experimental  curve  to  that  given  by  Eq.  (8),  one  car. 
attempt  to  "fit"  the  laboratory  curve  by  adjusting  K  and 
n  [14]. 

The  discussion  in  the  literature  centers  around  the 
fact  that  unique  K  and  n  cannot  be  found  to  fit  the  ex¬ 
perimental  curve  over  the  entire  range  of  oxygen  pressure 
from  zero  to,  say,  150  millimeters.  A  deterrent  to  fit¬ 
ting  the  curve  under  general  conditions  (i.e.,  physiologic) 
is  the  difficulty  of  determining  a  satisfactory  laboratory 
experimental  curve,  particularly  at  the  extreme  values  of 
oxygen  pressure  [14], 

But,  for  purpcses  of  this  comparison,  the  Dill  curves 
[16~!  are  satisfactory.  They  have  been  the  standard  for 
many  years  [7],  and  were  obtained  under  condition  reason¬ 
ably  well-controlled  for  the  time. 

Table  VIII  has,  in  the  red  cell  compartment,  a  reac¬ 
tion  of  oxygen  with  hemoglobin  (cf.  Table  V).  The  re¬ 
action  constant,  K,  is  determined  in  the  model  to  give  a 
normal  saturation  (97  percent)  at  100  mm  0^ .  Figure  l 
shows  two  cases  (computed  by  varying  the  CO2) :  two  curves 


Fig.  1— Experimental  hemoglobin  saturation  hypotheses 
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for  the  case  n  *  1  (as  in  Table  VIII)  and  a  curve  for  the 

case  n  =  2.5  as  Hill  originally  determined.  The  normal 

Dill  curve  also  appears.  The  Dill  curve,  of  course,  was 

determined  under  general  physiologic  conditions. 

Hemoglobin  has  four  binding  sites  for  oxygen.  So 

for  the  case  n  *  1,  in  order  to  bind  one  oxygen  per 

molecule  of  hemoglobin,  the  hemoglobin  was  conceptually 

broken  into  four  parts,  each  with  a  heme,  as  in  Eq.  (7). 

Thus,  for  the  case  n  *  1,  9.09  rather  than  2.2725  moles 

of  hemoglobin  combine  with  oxygen.  The  reaction  constant 

for  this  case  is  K  =  4.98x10  \  or  In  K  =  K'  =  -5.30. 

e 

Reducing  this  constant  immediately  decreases  the  maximum 

saturation  without  materially  shifting  the  curve  toward 

-3 

the  Dill  curve,  as  for  the  curve  with  K  =  6.74x10  , 

K'  =  5.0.  The  pH  is  7.39  in  plasma,  7.19  in  red  cells. 

The  Hill  curve  with  n  =  2.5  uses  the  proper  average 
number  of  moles  of  hemoglobin,  2.2725,  per  liter  of  blood, 
but  the  fixed  average  number  of  oxygen  molecules,  2.5, 
attach  to  each  molecule  of  hemoglobin.  Thus,  each  hemo¬ 
globin,  Hb^,  has  either  2.5  molecules  of  0^  attached  or 

none  at  all.  There  are  no  intermediate  stages.  An  equi- 

/ 

librium  constant  of  K  =  1.51x10  ,  K’  *  -11.1  is  necessary 

to  give  97  percent  saturation  at  100  mm  0^  and  the  same 


pH  as  before. 
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Some  of  the  data  for  Fig.  1  (the  n  *  2.5  and  n  -  1, 

K  ■  5.3  curves)  are  verified  in  Ref.  14  where  similar 
curves  were  computed  by  desk  calculator  (the  values  of 
K  are  not  given).  This  reasonably  validates  the  present 
model  under  varying  conditions.  But  the  discrepency 
between  the  empirical  Dill  curve  and  the  Hill  curves 
necessitates  better  hypothesis  and  an  improved  model. 
Subsequent  biochemical  hypotheses  for  the  oxygenation 
reactions  are  increasingly  sophisticated  [2]  and,  also, 
the  presence  of  CC^  and  H+  are  important.  Hemoglobin 
binds  CO2  in  the  carbamino  reactions  as  well  as  buffering 
the  H+  ion;  both  of  these  reactions  effect  oxygenation. 


Reference  2  discusses  these  problems. 
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Chapter  II 

CONCEPTUAL  COMPARTMENTS 


1 .  INTRODUCTION 


Chemically  speaking,  by  far  the  most  interesting 
constituents  of  the  blood  are  the  proteins.  These  species 
along  with  other  complex  const ituents--the  organic  phos¬ 
phates,  the  lipids,  the  carbohydrates--can  combine  with 
other  blood  constituents  potentially  to  create  innumerable 
biochemical  reaction  products.  However,  such  a  model,  with 
lists  of  possibly  thousands  of  reactions,  would  be  ex¬ 
tremely  unwieldy;  and  the  thousands  of  necessary  micro¬ 
scopic  reaction  constants  probably  could  not  be  determined. 

The  problem  then  is  to  make  a  model  which  approximates 
the  more  complex  real  situation  but,  at  the  same  time,  has 
exactly  the  same  characteristics  as  the  large  model  wherever 
possible . 

The  following  elementary  example  illustrates  this  re¬ 
mark.  Consider  a  polypeptide  having  two  identical  but  in- 
depent’ei  t  ionizing  groups- - i . e . ,  the  ionization  of  one  does 
not  affect  the  ionization  of  the  other.  Thus,  a  representa¬ 
tive  molecule  HRH  has  three  ionization  states  and  four 


reactions 
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HRH  -  HR-  +  H+  ,  kL 
HRH  -  -RH  +  H+  ,  k2 

(1) 

HR-  -  -R-  +  H+  ,  k3 
-RH  -  -R-  +  H+  ,  k4 

each  with  the  intrinsic  microscopic  dissociation  constant 
k.  o’:  incrinsic  association  constant  1/k,  .  But  if  the  paths 
of  ionization  are  inmaterial,  one  ran  instead  write  [8] 

HRH  -  (HR-  +  -RH)  +  H+  ,  ^ 

(2) 

(— RH  +  HR-)  -  (— R— )  +  H+  ,  K2 

with  constants  and  K^,  respectively,  where  and 
are  the  apparent  titration  constants  obtained  from  a 
Laboratory  titration  curve,  and  HR—  and  — RH  are  indis¬ 
tinguishable,  i.e.,  a  single  species.  The  constants  k^ 
and  can  be  related  by  the  following  reasoning:  Con¬ 
sidering  the  several  mass  action  equations  for  Eqs.  (1) 
and  (2),  we  may  write,  as  in  Ref.  8, 
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v  [(-RH)  ±  tHR-)]  (H+) 
1  (HRH) 


HUQ  (H+)  ,  (HR-)(H+) 
(HRH)  (HRH) 


(3) 


+  k 


2 


'1nd 


.  _ (-»-?  (»* 

2  [(-RH)  +  (HR-)' 


)  [(-RH)  +  (RH— )  1  i 
(— R—  (H+)  [ 

^  i 


(4) 


The  two  chemical  systems,  Eq.  (1>  and  Eq.  (2),  are  thus 
completely  convertible,  and  equivalent  in  the  sense  that 
the  computed  ratios  of  concentrations  for  (— R— )/(HRH)  at 
equilibrium  are  identical,  as  are  the  concentrations  of 
(H+)  and  (-HR  +  HR-).  The  first  case,  however,  requires 
four  reactions  to  represent  the  system;  the  second,  only 
two . 

Tills  and  other  types  of  reductions  to  be  discussed 
have  a  commc,s  'perty:  They  reduce  the  s  ize  of  the 
problem  but  not  the  simulation  accuracy.  Such  reductions 
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are  important  In  the  computer  model,  since  both  the  com¬ 
puting  time  and  the  size  of  the  computer  required  for 
a  solution  increase  nonlinearly  with  the  problem  size. 

In  biological  systems,  a  problem  could  easily  be  too 
large--!. e. ,  have  too  many  chemical  reactions  in  its 
representation-- for  even  a  reasonably  large  computer. 

In  the  blood,  serum  albumin  has  190  sites  for  active 
hydrogen  ionization,  each  of  which  falls  into  one  of  four 
classes.  These  classes  are  independent  but  not  mutually 
exclusive;  i  c. .  ionization  in  one  class  does  not  prohibit 
ionization  ir.  another.  Because  ionization  can  occur 
simultaneously  in  all  four  classes,  they  must  be  con¬ 
sidered  as  simultaneous  events  though  not  with  equal 
probability.  Since  each  site  has  two  possible  states  and 
there  are  99  sites  in  the  first  class,  16  in  the  second, 

UQ  \  f; 

57  in  the  third,  and  13  Ln  the  fourth,  then  2  "  *  2  » 

2^7  y  2^  »  1060  states  are  possible--each  of  which  must 
be  represented  in  the  model.  Obviously,  some  reduction 
Is  necessary. 

While  an  equivalent  representation  can  reduce  the 
size  of  complex  systems,  as  in  £qs  .  (1)  and  (2).  con¬ 
ceptual  eompartroentai izat ion  is  generally  an  even  more 
economical  device.  As  will  be  shown,  the  conceptual 
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compartments  are  artificial  compartments  created  in  the 
chemical  milieu  of  the  model  to  isolate  classes  of  re¬ 
actions  on  complex  molecules  when  these  reactions  can  occur 
simultaneously.  As  in  the  example  of  serum  albumin,  the 
reactions  of  independent  classes  may  occur  simultaneously-- 
thus,  the  same  molecule  must  have  multiple  representation 
to  show  all  possible  configurations.  The  procedure  sug¬ 
gested  here  is  to  identify  and  characterize  the  myriad 
chemic.  <  reactions  according  to  class  (type  of  ligand) 
and  interdependence  (mutually  exclusive,  dependent,  or 
independent),  in  ;>rder  to  represent  each  c lass  of  reaction 
as  an  isolated  event  rather  than  representing  all  possible 
td-class  species.  Since  the  classes  of  reactions  are 
much  fever  than  the  total  possible  number  of  species,  the 


reduction  can  be  considerable. 
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2.  CONCEPTUAL  COMPARTMENTS 

PROTEIN  IONIZATION- -EXAMPLE 

The  blood  macroroolecules  of  principal  concern  are 
serum  albumin,  the  various  globulins,  and  hemoglobin. 
Because  serum  albumin  is  perhaps  the  simplest,  most 
stable,  and  best  known,  it  will  be  used  extensively  as 
an  example.  This  protein  is  the  most  abundant  of  the 
plasma  proteins,  about  ten  times  more  numerous  than  the 
next  most  abundant.  The  pertinent  characteristic  of 
serum  albumin  is  its  ability  to  buffer  hydrogen  ion  in 
the  physiological  range  of  pH,  but  its  activity  is  by  no 
means  limited  to  proton  binding.  Serum  albumin  also  binds 
small  ions  of  either  sign  and  even  small  neutral  molecules. 
To  simplify  the  description  of  serum  albumin  for  the  model, 
we  will  deal  exclusively  with  thw  hydrogen  binding  activity. 

Tanford,  et  al,  [17]  and  Tanford  [18]  (sources  for 
much  of  the  factual  data  on.  albumin)  list  among  the 
hydrogen-binding  sites  of  serum  albumin.  99  carboxyl 
radicals  with  average  pK  •*  4.02  at  25°C.  The  carboxyl 
groups  alone  have  99  configurations  of  the  eolecule  in 
which  it  hat  lost  just  one  hydrogen,  99*98/2  configura¬ 
tions  in  which  it  has  lost  two  tydrogens .  and  so  on,  in 
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99 

all  2  species  in  Che  milieu.  But  several  ways  exist 
Co  reduce  this  variety  to  an  approximate  or  an  equivalent 
system. 

Assume  first  that  the  99  sites  are  equivalent  and 
independent,  that  the  ionization  of  one  site  will  not 
affect  the  probability  of  ionization  of  any  other.  An 
equivalence  can  then  be  drawn  between  the  ionization  of 
this  polybasic  acid  and  99  monobasic  acids,  each  with  the 
same  pK  and  the  same  concentration  as  the  protein.  But 
this  equivalence  is  immediate  if  the  protein  is  concep¬ 
tually  split  into  n  parts.  If  h  is  the  number  of  moles 
of  hydrogen  released  per  mole  of  protein  per  site  then, 
as  for  a  monobasic  acid, 

RH  -  R-  +  H+  ,  k  (5) 


or 


(RH) 


(H+) 


1-h 


<e>) 


where  k  la  the  microscopic  dissociation  constant  for  that 
particular  site.  If  this  were  the  only  site  for  ioniza¬ 
tion,  then  clearly  the  ideal  titration  constant  determined 
in  the  laboratory  would  be  k.  For  n  identical,  independent 
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sites  on  each  molecule  of  R,  n  times  as  many  protons  will 
be  released  at  the  same  pH,  but  of  course  the  ionization 
pK  remains  the  same,  l.e.,  -log  k.  Also,  u  *  nf\  is  the 
moles  of  protons  for  all  sites  per  mole  of  R.  Therefore, 


nil 

n(l-h) 


IB-1 

(RH) 


or 


n-  v 


where,  again,  the  empirically  observed  titration  constant 
in  vitro  would  also  be  simply  k. 

For  the  similar  case  of  n  moles  of  monovalent  acid, 
each  with  the  same  concentration  and  pK  as  for  R,  the 
titration  constant  will  equal  that  of  the  protein,  i.e., 
either  system  will  buffer  maximally  at  the  pK  of  the  in¬ 
trinsic  ionization  groups.  Then,  in  titration  or  buffering 
action.  Che  two  systems  would  be  indistinguishable  except 
for  electrostatic  concentration  effects. 
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We  note  for  future  reference  the  difference  between 
the  titration  constants  and  the  dissociation  constants. 

For  the  monovalent  acids,  these  two  constants  are  the 
same  and  are  given  by  the  mass  action  equation,  l.e.,  the 
pK  of  the  ionizing  hydrogen.  For  the  protein,  the  pK  of 
the  ionizing  class  (Eq.  (6))  gives  the  titration  constant; 
but  the  dissociation  constant,  given  by  mass  action 
equations  where  the  protein  is  not  conceptually  split, 
depend  upon  the  number  of  hydrogens  dissociated  (and 
perhaps  the  particular  configuration  of  their  dissocia¬ 
tion).  This  difference  is  important  because  the  apparent 
mass  action  equation  constant  for  a  protein  may  be  dif¬ 
ficult  to  measure,  whereas  the  titration  constant  can  be 

♦ 

determined  from  the  in  vitro  titration  curve. 

This  discussion  implies  that,  under  the  strong  assump¬ 
tion  of  equivalence  and  independence  of  the  groups  of  e 
piven  class,  for  application  to  the  computer  model,  the 
simulation  of  the  protein  could  be  obtained  by  writing  down 
the  chemical  reaction,  Eq.  (3).  and  inserting  this  in  the 
model  with  constant  k  and  in  amount  n  times  RH,  i.e.,  n 
monovalent  acids.  In  fact,  the  simulation  is  not  quite 
so  simple:  the  milieu  would  contain  n  times  too  much 
protein  and  the  osmolarity  would  be  in  error. 

*Ref.  8,  Chap,  9. 
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However,  in  other  re8pects--e.g. ,  the  buffering  re- 
«ction--this  simulation  would  be  perfectly  valid.  There¬ 
fore,  the  model  may  reasonably  contain  an  osmotically 
isolated  compartment  for  the  monovalent  acids.  This 
conceptual  compartment  will  be  impermeable  to  all  species 
except  H+  ion  and  will  not  be  counted  in  the  osmolar ity-- 
it  will  be  a  sink  or  source  of  H  ion,  a  place  to  store 
H+  Ion  effectively  out  of  the  milieu  in  proportion  to  the 
pK  of  the  monovalent  acids  which  are  availaole  only  with¬ 
in  the  conceptual  compartment. 

The  example  of  serum  albumin  in  Sec.  1  permits  more 
specificity.  Table  X  gives  five  major  classes  of  ionizing 
sites  on  the  protein,  each  with  a  distinct  number  of  sites. 
There  are  17  imidazole  groups,  57  c-amino  groups,  19  phenolic 
groups,  and  22  guanidine  groups,  so  that,  as  on  p.  46, 
approximately  10^  species  are  possible  and  an  equal  number 
of  intrinsic  constants  could  be  required  for  a  chemical 
description. 

The  method  of  conceptual  compartments  replaces  each 
ionizing  class  with  a  single  monobasic  acid  having  the 
"average"  pK  of  that  class,  but  with  n  times  the  mole 
number  of  the  protein,  where  n  Is  the  number  of  sites  in 
the  class.  Then,  to  not  disturb  the  osmolar ity  (or  the 
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Table  xa 

TITRATABLE  GROUPS  PER  65,000  GRAMS  SERUM  ALBUMIN 


Found  by 
Titration 

Amino  Acid 
Analysis 

orcarboxyl 

(1) 

1 

0,  y-carboxyl 

99 

*89.9 

Imidazole 

16 

16.8 

At- ami  no 

(1) 

1 

c- amino 

57 

57.0 

Phenolic 

19 

18.1 

Guanidine 

22 

22.0 

Sulfhydryl  (free) 

(0) 

0.7 

3 

Adapted  from  Tanford,  et  al .  [17". 


total  mole  count),  an  artificial  compartment  is  created 
to  contain  the  monobasic  acid.  This  compartment  has  but 
two  species,  the  ionized  and  un- ionized  forms.  The  dis¬ 
tribution  between  these  two  species  is  determined  by  the 
pK  of  the  acid  and  the  ion  concentration  of  the  liquid 
compartment . 

Thus,  the  serum  albumin  problem  is  reduced  to  a 
representation  containing  just  eiever  species:  the  protein 
in  the  main  compartment  and  two  species  in  each  ionization 
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class  compartment.  The  remainder  of  this  section  shows 
in  what  respect  this  problem  with  eleven  species  (plus 
the  remaining  species  in  the  liquid  phase,  such  as  Na+, 

Cl  ,  OH  ,  ^0,  i  j.)  is  equivalent  to  the  original  problem 
of  10^  species,  where  "equivalent0  means  that  the  solu¬ 
tions  of  these  two  problems  will  be  indistinguishable. 

Shapiro  <  4  1  develops  and  proves  in  general  an 
equivalence  theorem  adaptable  to  the  particular  case  of 
serum  albumin.  Suppose,  first,  that  the  ionizing  sites 
on  the  protein  are  divisible  into  classes  (e.g.,  carboxyl, 
alpha-amino,  etc.)  and  that  within  each  class  each  ionizing 
site  is  indistinguishable,  the  order  of  ionization  im¬ 
material,  and  the  free  energy  associated  with  ionization 
is  strictly  additive  (i.e.,  the  ionization  of  one  site 
does  not  in  any  way  affect  the  Ionization  of  any  other). 

Second,  suppose  this  entire  protein-carboxyl  system 
is  replaced  in  the  mathematical  model  with  a  single  mono¬ 
basic  acid  having  a  pK  identical  to  the  equivalent  hydrogens 
of  the  protein  but  99  times  the  number  of  moles  of  the 
protein.  Preserving  the  osmolality  requires  a  new  phase 
or  conceptual  compartment  containing  only  the  monobasic 
acid  in  its  two  states,  ionized  and  un-lonized  the  pro¬ 
tein  itself  remains  unreactive  in  the  original  milieu 
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(again,  to  keep  the  osmolality  correct).  The  protein  and 
the  monobasic  acid  are  related  by  a  conservation  of  mass 
equation. 

Shapiro  then  shows  that  under  these  conditions  the 

total  number  of  moles  of  ionized  monobasic  acid  will 

equal  the  number  of  ionized  sites  in  the  original,  large 

problem,  although  the  number  of  species  has  decreased 
99 

from  2  to  just  two.  The  retraining  species  other  than 
the  protein  and  monobasic  acid  will  be  unaffected  in 
concentration  or  distribution.  Further,  the  total  thermo 
dynamic  free  energy  of  the  system  will  be  unaffected.  We 
forego  the  details  or  proof  of  the  more  general  theorem. 

This  model  then  proceeds  similarly  for  the  other 
classes  of  ionising  sites  on  the  protein  (Table  X).  The 
pK  of  each  class  comes  directly  from  Tanford,  et  al .  [171 
who  determine  the  pK  from  the  titration  curve  with  a 
correction  for  the  static  charge  on  the  molecule  and  for 
the  adsorbed  chloride  ion.  Reference  19  discusses  this 
computer  model  in  detail  and  shews  thfit  it  reproduces 
the  titration  curve  for  scrum  albumin  within  two  ionizing 
sites  over  the  range  of  pH  from  3.5  to  10.  The  extreme 
ends  of  the  curve  require  a  more  complex  model. 
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3.  BLOOD  PROTEINS:  PLASMA 

The  power  and  flexibility  of  conceptual  compartmen* 
talization  permits  further  development  of  the  elementary 
and  insufficiently  accurate  blood  model  of  Chap.  I.  The 
simulation  will  continue  by  incorporating  further  com” 
plexities  of  the  hemoglobin  reactions  and  certain  bu£~ 
fering  systems  of  the  blood. 

■*tn  one  view,  the  principal  buffering  system  of  the 
blood  the  bicarbonate.  With  the  loss  of  each  CC^ 
flTv’ecuie  at  the  lung  surface,  a  free  hydrogen  ion  is  ab¬ 
sorbed  in  solution  by  the  OK  ion  from  the  decomposing 
bicarbonate  ion,  thus  making  the  blood  more  basic  At 
the  tissue  level,  the  absorption  by  the  blood  of  each 
molecule  creates  hydrogen  ions,  making  the  blood  more 
ac  id . 

These  two  systems-* the  production  and  consumption  of 
H  ion  by  the  bicarbonate  react  ions** need  not  be  and  per* 
haps  seldom  are  in  ex'1***  balance,  thus  causing  respiratory 
acidosis  or  alkalosis.  However,  the  vrery  low  CQ0  content 

m, 

of  the  atmosphere  relative  to  the  blood  permits  CO,  to  be 

£ 

easily  <c  leased  *2,  tne  Kings;  the  autonomic  controls  over 
deptn  and  race  cf  ^ esnii at  ion  and  blood  flow  rate  tend  to 
balance  production  ard  absorption  of  i^n. 


The  necessity  for  this  buffering  may  be  eisily  shown. 
At  a  steady  resting  state,  blood  coursing  through  the 
tissues  picks  up  about  40  ml  of  CO2  Per  liter  in  changing 
from  arterial  to  venous.  This  is  equivalent  to  the  re* 
lease  of  about  2  mM  of  hydrogen  ion,  since  C0o  either 
combines  with  ^0  to  form  HCO^  and  liberate  a  H+  ion, 
or  forms  a  carbamino  compound  with  hemoglobin  and  again 
liberates  almost  one  hydrogen  per  reaction.  Since  the 
hydrogen  ion  concentration  is  normally  only  about  10 
moles  per  liter  of  blood,  this  tremendous  excess  acidity 
must  be  buffered.  The  net  change  in  pH  during  the  normal 
transformation  from  arterial  to  venous  blood  is  only  about 
0.03  pH  units. 

Hemoglobin  itself  takes  up  approximately  three- 
quarters  of  the  excess  hydrogen  in  the  oxylabile  binding 
of  hydrogen  on  group*  which  become  available  on  the  molecul 
as  the  oxygen  is  removed  at  the  tissue  level.  The  oxygen¬ 
ation  of  hemoglobin  yields  this  hydrogen  at  the  lung  level, 
thus  preventing  the  blood  from  going  basic  from  the  loss 
of  C0r 

The  remaining  excess  hydrogen  at  the  tissue  level  is 
buffered  principal iy  by  the  plasma  proteins,  the  phosphate 
ionization  systems,  and  by  the  non-oxvlabiie  sit ~ on  the 
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hemoglobin  [20].  Plasma  proteins,  however,  may  account 
for  a  maj ^r  portion  of  the  remaining  excess,  since  the 
imidazole  groups  of  the  serum  albumin  are  at  their  maxi¬ 
mum  buffer  inf  power  in  this  important  range  of  pH  of  the 
viable  system,  referring  again  to  Tanford  [17],  one 
millimole  of  serum  albumin  buffers  approximately  one 
millimole  of  hydrogen  in  a  change  of  0.1  pH  point.  Since 
i;he  pH  changes  about  0.03  points  from  arterial  to  venous 
blood  and  since  there  is  approximately  0.5  mM  of  serum 
albumin  per  liter  of  blood,  the  serum  albumin  could  buffer 
up  to  0.15  mM  of  H+  during  this  change. 

incorporating  the  plasma  protein  buffering  system 
into  the  mathematical  model  requires  modifying  the  gen¬ 
eralized  model  of  plasma  protein  described  in  Chap.  I. 
Serum  albumin  accounts  normally  for  80  or  90  percent  by 
mole  number  of  the  total  plasma  protein,  nd  the  various 
globulins  for  most  of  the  remainder.  Fibrinogen  is  also 
present,  but  i>  n^t  thought  to  be  an  importer  buffer. 

We  assumed  that  the  two  principal  proteins  represent  all 
Qt  the  important  buffering  power  of  the  plasma  fraction, 
and  have  similar  titration  curves,  particularly  in  the 
middle,  viable  range  o(  pH.  This  latte**  assumption  is 
not  quite  true  '18,21],  but  considering  the  small 
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proportion  of  globulin  and  the  reasonably  small  dif¬ 
ference  in  the  titration  curves  in  the  viable  range,  the 
error  will  be  insignif icant--not  more  than  2  percent  of 
the  maximum  buffering  power  of  the  plasma  proteins. 

Finally,  the  titration  curve  itself  is  more  sophis¬ 
ticated  than  indicated  in  the  previous  section.  Tanford 
[17]  used  a  solution  of  serum  albumin  with  salt  added  to 
provide  an  ionic  strength  0.15.  Thus,  the  titration  curve 
incorporates  the  electrostatic  charge  on  the  molecule  of 
protein  and  the  comne  tive  binding  of  chloride  ion.  The 
pK  of  the  hydrcgen-binding  sites  thereby  varies  appro¬ 
priately  through  the  range  of  pH.  References  19  and  22 
detail  this  problem.  Of  course,  the  behavior  of  the 
protein  in  plasma  may  differ  greatly  from  its  behavior 
in  the  0.13  molar  NaCl  titrating  solution.  However,  this 
error  can  temporarily  be  discounted,  since  Chap.  Ill  will 
demonstrate  with  Astrup  diagrams  and  more  stringent  tests 
Chat  the  simulated  biood  does  buffer  properly. 

Serum  albumin  binds  anions,  cations,  and  even  un¬ 
charged  particles  by  incompletely  understood  mechanisms  '8 
However,  Cl  probably  binds  to  the  positively  charged 
amine  groups;  this  reaction  is  a  function  of  the  pH  and 
the  concentration  of  chtoride  ions.  Scatchard.  et  al .  [23 
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for  example,  discuss  the  binding  of  Cl  ,  which  is  im¬ 
portant  in  the  respiratory  function  of  blood;  Loken, 
et  al,  [241  treat  the  binding  of  calcium.  These  reac¬ 
tions  change  the  net  charge  on  the  protein  molecule  and 
thus  affect  the  buffering  power  of  the  protein;  but  by 
tying  up  permeable  ions,  they  also  affect  the  Gibbs- 
Donnan  equilibrium.  Neither  of  these  systems  are  ex¬ 
plicitly  in  effect  in  the  model,  where  serum  protein  may 
bind  up  to  50  percent  of  the  calcium,  and  up  to  10  ions 

M* 

of  Cl  per  molecule  of  protein,  depending  upon  the  pH. 

The  binding  of  these  Lons  by  other  proteins  has  not 
been  studied  as  closely;  in  particular,  hemoglobin's  role 
is  net  well  understood ,  We  presume  that  hemoglobin  and 
other  complexes  vithir  the  cell  may  also  bind  these  small 
ions.  Although  we  may  assume  that  none  of  the  single- 
valent  alkaline  cations  are  bound,  both  the  double-valent 

■4-f 

Hg  and  the  double-valent  anions  are  tied  up  in  various 
ways.  The  literature  regarding  these  intracellular  re¬ 
actions  is  not  as  definite  as  for  Cl”  and  Cs*^,  and  no 
attempt  has  been  made  to  account  Cor  them. 


i 


H 
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4.  BLOOD  PROTEINS:  RED  CELL 


Hemoglobin  Is  only  one  Important  constituent  of  the 
red  cell.  In  fact,  the  water  content  is  approximately 
99.45  percent  of  the  total  moles;  hemoglobin  is  only 
about  1.5  percent  of  the  remainder .  or  about  2.30  mM  per 
liter  of  blood.  However,  hemoglobin  is  more  prominent  in 
terms  of  ma3s:  each  liter  of  blood  contains  about  160 
grams  of  hemoglobin- -all  in  the  450  cc  occupied  by  red 
cells.  The  hemoglobin  protein,  minor  impermeable  con¬ 
stituents,  and  the  membranes  of  the  red  cells  (which  are 
very  chin)  increase  the  volume  of  the  cells  over  that  of 
the  pure  water  content  by  about  120  cc  per  liter  of  blood. 

The  physical-chemical  structure  of  the  red  cell 
Interior  may  orient  the  hemoglobin  molecules  with  respect 
to  each  other,  perhaps  linked  in  a  fluid  crystalline 
lattice.  In  any  case,  we  assume  an  homogeneous  fluid 
throughout,  with  reasonably  steady  composition  states 
(within  measurable  error)  when  the  cell  environment  is 
constant.  A  change  in  the  environment  results  in  a  change 
in  the  cell  interior — l.e..  In  the  rate  of  it*  meci  misres-- 
end  consequently  a  new  steady  state. 

Besides  hemoglobin,  inorganic  electrolytes,  organic 
permeable  molecules  and  ions,  and  water  the  red  cell 
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contains  a  complex  of  interesting  and  active  substances, 
most  of  which  are  charged.  The  kinds  of  these  particles 
(organic  phosphates,  enzymes,  chelates,  etc.)  are  so 
numerous  and  their  specific  contributions  so  varied  that 
they  are  usually  lumped  together  in  the  electrolyte  balance 
as  simply  X  . 

The  milliequivalents  of  this  total  are  more  easily 
calculated  than  the  actual  millimoles.  By  assuming  that 
the  red  cell  milieu  of  a  given  volume  of  arterial  red 
cells  is  neutrally  charged,  one  can  compute  the  total 
charge  attributable  to  X  --about  10  milliequivalents  of 
negative  charge  per  liter  of  blood.  Then,  since  the 
osmolarity  of  the  milieu  should  equal  that  of  plasma, 
and  the  moles  of  each  other  substance  are  report _d,  one 
can  determine  that  there  should  be  about  7.5  millimoles 
of  X,  only  an  approximate  calculation. 

From  this,  however,  the  average  charge  per  average 
molecule  is  between  one  and  two.  about  the  average  charge 
of  inorganic  phosphate  or  of  a  hemoglobin  molecule  near 
the  viable  pH.  Accordingly,  the  mathematical  model  carries 
7.6  raM  of  X  per  liter,  but  (as  Chap.  Ill  proves)  for  the 
present  model  the  net  charge  need  not  be  specified.  This 
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is  a  temporary  procedure  for  handling  the  mixture  of 
substan  is  denoted  by  X  ,  whose  details  may  be  added 
gradually  to  the  model  as  specific  reactions  data  become 
available  or  important.  However,,  the  situation  for 
hemoglobin  differs  considerably,  since  its  details  are 
important  and  relatively  well-known. 

Hemoglobin's  chemical  reactions  important  to  respira¬ 
tion  are  the  reversible  combinations  formed  with: 

1)  Oxygen  (via  the  four  i'.ron  atoms  in  the  heme 
groups); 

2)  Carbon  dioxide  (via  ~Nh ,,  groups,  probably  those 
in  heme  and  lying  near  the  iron  atom) ; 

3)  Hydrogen  ion  (via  the  -COOH,  and 

groups,  a  significant  number  of  which  are  oxy- 
lanile  in  their  ionization;  i.e.,  the  oxygen 
attachment  affects  the  bonding  energy). 

The  Adair  [25 1  equation  for  the  oxygen  dissociation 
of  hemoglobin  under  standard  conditions--viz , , 


,  (atp  +  2a,p2  ♦  3a3p3  f  4a4pA) 

ico  "  77  r  I  H  T~  V 

4(aQ  ♦  ♦  a2p  ♦  a^p  ♦  a^p  ) 


(7) 


where 
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y'  ■  percent  saturation  of  hemoglobin, 


ao  ■  ^ 


TTk. 


p  -  oxygen  pressure  in  millimeters, 


CHb^O.,)^  -  Kt  PrHb4(02)l.1 


--is  used  in  this  model,  as  in  Ref.  2,  for  the  generation 
of  the  empirical  saturation  curve. 

Given  an  empirical  curve,  y  vs .  p,  from  the  labora¬ 
tory,  the  analytic  problem  is  to  find  the  unique  set  of 
four  parameters,  the  ratios  a  :  a.:  a^:  a  •  a.,  or  simply 
the  apparent  ,  that  enable  Eq.  (7)  to  reproduce  the 
curv-~  within  the  model.  Because  of  the  oxylability  of 
the  ancillary  reactions  of  hemoglobin,  however,  these 
apparent  may  differ  from  the  intrinsic  of  Eq .  (9). 

In  this  circumstance,  Pauling  [26^  correctly  wrote 
for  the  apparent  constants  in  the  milieu  and  for  the 
intrinsic  constants  of  “pure”  hemoglobin.  The  solution 
of  this  impasse  lies  in  finding  analytical  transformat  ions , 
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( 


based  upon  structural,  physical,  and  chemical  hypotheses, 
which  will  carry  one  set  of  constants  into  ‘he  other, 


kj  "Kr 


The  apparent  macroscopic  KJ , 


HM02>j-i  +  °2  r  HWj  £or  each  j  ’  <10> 


are  not  identical,  and  we  have 


-lil 

100 


K^p  +  2K[K^p2  +  3KjK£K^p3  +  4K[K^K^p4 
4  +  K* p  +  KjK^p2  +  K-jK^K^p3  +  KjK^Klp4 


(ID 


12  3* 


1  2  3  4* 


for  the  Adair  Eq.  (7). 

Reference  2  discusses  the  use  of  the  computer  in  the 

derivation  of  the  oxygenation  constants.  This  study 

utilizes  the  ratios  derived  independently  on  the  computer 

and  by  Roughton,  at  al.  [27,28',  namely,  Kj,  :  *2  *  *3  * 

K!  : :  l  :  0.4  :  0.24  :  9.4.  The  absolute  values  are  «d- 
4 

Justed  to  give  the  proper  saturation  of  hemoglobin  at  p50, 
the  oxygen  pressure  for  half-saturation.  Chapter  III 


-66- 


discusses  the  quantitative  goodness-of-f it  for  these 
coni?  tants . 

This  simulation  assumes  that  the  transformation  T 
is  affected  only  by  the  ancillary  oxy labile  carbamino 
and  hydrogen  ionization  reactions.  In  deriving  T,  we  may 
reasonably  assume  [141  two  oxy labile  hydrogen  ionization 
sites  per  heme  group  on  hemoglobin,  and  16  sites  per  heme 
for  non-oxylabile  ionization.  Identification  of  the 
position  and  kind  of  these  sites  on  the  actual  molecule 
is  not  yet  definite  [14],  but  for  convenient  reference 
this  simulation  classifies  the  oxylabile  sites  as  a 
carboxyl  and  an  alpha-amino,  the  non-oxylabile  as  iden¬ 
tical  ligands  on  histidine  groups.  Further,  we  assume 
that  the  pK's  of  the  oxylabile  ^roup3  are  sufficiently 
far  apart  that  only  three  species  with  respect  to  the 
oxylabile  heme  groups  appear  in  the  milieu: 


R 


COOH 


K 

J  coo" 

R 

nh3 


K. 


COO 


(12) 


where  R  is  a  heme  group. 

Evidently,  the  CO^  molecule  is  bound  only  by  the 
amphat.  he  spec  is?  on  the  rf«ht  of  Eq .  il2),  *rfith  the 
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amino  in  the  NH2  condition.  The  apparent  constants  for 
oxygenated  and  reduced  hemoglobin  differ,  and  therefore 
the  equations  for  the  model  are 


Hb  +  C02  ~  Hb(C02)  +  H+  ,  Krc  -  6,4xl0*6  (13) 

and 

Hb02  +  C02  -  Hb02(C02)  +  H*  ,  KQ(,  -  2.5xl0‘6  (14) 

where  the  symbol  Hb  Indicates  a  1/4  hemoglobin  molecule, 
and/or  a  reactive  group  on  an  amphan ionic  heme.  Sim¬ 
ilarly,  the  pK^  for  hydrogen  ionization  will  have  sub¬ 
scripts  referring  to  oxygenated  or  reduced  hemoglobin, 
and  the  apparent  pK^  may  correspond  to  any  number  of  H+ 
equivalents  ionized.  The  equations  relating  to  the  first 
stage  of  oxygenation  are: 


Hb,  +0,  -  Hb,  0-  .  K! 

4  i  4  2  1 

Hb*  r  Hbt  ♦  H*  •  *  b  2b 

HbJ  *■  Hb4  ♦  *  P\i  *  7.93 


Hb4°2  r 

Hb^O*  4  H*  , 

<>K0Z 

-  5.73 

uk  n*  r 

»  U * 

**  2 

♦  ♦ 
uk  n  ±  ii 

»  6 .  hfi 

(15) 
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where  In  the  first  equation  is  the  int7:inoic  reaction 

V 

constant  for  the  attachment  of  0^  on  any  of  the  four  hemes. 
KJ  is  the  apparent  constant  measured  in  the  laboratory  in 
the  presence  of  H*  or  C(>2  reactions.  Equations  (13),  (14), 
and  (15)  specify  the  Interrelated  reactions  with  respect 
to  the  first  oxygenation  state.  Computing  the  ratio  (in 
moles)  of  oxygenated  to  reduced  molecules  yields  (sup¬ 
pressing  the  square  concentration  brackets) 


Hb,  0.. 
2 

Hb. 

4 


W 


Hb402  +  Hb402  >  Hb40‘  +  Hb^CCQ,) 
Hb*  +  Hb*  +  Hb"  +  Hb4(CO„) 


(16) 


Alternatively,  using  the  equilibrium  constants  from  the 
individual  mass  action  equations,  Eqs.  (9).  (13),  (14), 
and  (15),  pOj  in  millimeters,  and  concentrations  of  each 
species  (suppressing  the  square  brackets),  yields 


Hb  0, 

~Hb~  '  V 

4 


4  V  t* 

H  „  01  K 


K  p  +  K  n  - _  -f  K  n  il  *  >  t-  2  01 

l  ’  r  IL.,  lp  ♦  r  ~  K 

_ H  1  H  H 


or. 


I  ♦ 


H 


H 


St  I  ^  C°2  K0I 

“T7  *c 
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where  numerator  and  denominator  are  divided  by  fHb,  "!  and 

4 

multiplied  by  760  mm.  Thus,  a  complex  ratio  modifies  the 
intrinsic  constant  to  give  K',  the  apparent  constant 
for  first-stage  oxygenation  in  the  presence  of  ancillary 
reactions.  This  result,  Eq.  (17),  resembles  derivations 
elsewhere,  e.g.,  Pauling  '26j. 

Equation  (11),  corresponding  to  the  Adair  formulation 
for  the  dissociation  of  hemoglobin--!. e. , 


ajp  4  Sa^p*  *  3a^p3  ♦ 

♦  »{P  ♦  *2p2  *  a3p3  * 


(18) 


shows  the  constant  a’  -  K’.  Now,  if 


7 
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a'2  >  K*  , 

a3  *  K[  K2  K3  ’ 

a!  -  K’  Ki  Ki  K!  , 

4  12  3  4 

tH 

where  Che  K'^  are  the  apparent  j  -stage  oxygenation  con¬ 
stants  in  the  milieu  of  reactions,  then 


an  argument  similar  to  that  yielding  Eq.  (17). 


This  mathematical  model  incorporates  Eqs.  (19).  A 
different  model  may  result,  hovevtsr,  if  the  relations  of 
the  several  reactions  of  h««oglobin  are  eventually  <je- 

4* 

ter mined  to  be  dependent  in  other  ways.  Some  of  the  H 
ionization  must  be  oxylabile,  but  certainly  not  all;  and 
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the  oxygen  absorption  will  depend  on  H+  lability.  Simi¬ 
larly,  remote  areas  of  the  protein  and  heme  groups  may 
bind  CO2  as  non-oxylabile  carbamino,  although  carbamino 
formation  always  depends  or.  the  pH.  Also,  some  of  these 
reactions  may  be  mutually  exclusive--e.g. ,  the  CO  and  0^ 
reactions,  and  the  Ca  and  H  reactions,  not  considered 
here. 

In  the  format  of  the  mathematical  model,  each  stoi¬ 
chiometric  chemical  reaction  requires  a  separate  entry  in 
the  matrix  of  reactions.  Thus,  if  five  ways  exist  for 
O2  to  combine  with  Hb^  and  four  ways  for  H+  ionization, 
and  if  these  reactions  are  not  mutually  exclusive,  then 
20  entries  would  represent  the  possible  combinations. 
However,  conceptual  compartments  convert  this  "multipli¬ 
cative"  model  into  an  equivalent  "additive"  model  of  nine 
entries.  In  addition  to  the  principal  red  cell  compart¬ 
ment,  the  latter  model  computes  carbamino  reactions  and 
the  H*  ionization  of  reduced  and  oxygenated  hem*  in 
separate  but  dependent  conceptual  compartments. 

Equations  (19)  relate  to  a  model  which  compute?  the 
carbamino  reactions  of  oxygenated  or  reduced  heme  the 
corresponding  ionization  compartments .  This  procedure 
is  efficient  and  does  not  yield  multiplicative  combinations 


because  Hb+,  Hb*,  Hb” ,  and  HbCC^  are  assumed  co  be 
mutually  exclusive  species,  CO2  combining  only  with  the 
amphanion  form.  Also,  the  present  model  regards  the 
product  of  Eqs.  (13)  and  (14)  as  completely  ionized,  an 
assumption  which  is  probably  not  much  in  error  at  viable 
pH. 

Finally,  the  computer  program  incorporates  Eqs.  (9), 
(13),  (14),  and  (15)  as  lists  of  the  expected  chemical 
reactions  and  species  in  appropriate  compartments.  Table 
XI  shows  the  red  cell  and  subsidiary  conceptual  compart¬ 
ments  for  the  current  mathematical  model.  In  the  red 


cell ,  the  Adai’'  formulation  has  four  stages  of  oxygena¬ 
tion  of  hemoglobin.  The  coefficients  of  the  terms  re¬ 
duced  heme,  oxy  heme,  and  histidine  give  the  number  of 


equivalents  assumed  in  n 


j 


moles  of  hemoglobin  for  sub¬ 


sequent  reactions  in  the  conceptual  compartments.  Thus, 


the  coefficients  in  the  red  cell  compartment  count  the 


moles  of  reduced  and  oxygenated  hemes  and  the  number  of 
oxystable  H  ionization  sites,  which  are  assumed  to  be  on 


histidines.  The  two  subsidiary  cc  apartments  each  contain 
the  two  stages  of  Ionization  and  a  carbamino  reaction. 

The  amounts  of  reduced  heme,  oxy  heme,  and  histidine 
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Extracted  from  Tablr  XIII,  p.  102. 
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species  in  the  conceptual  compartments  derive  from  the 
red  cell  via  conservation-of-raass  equations  called 
"transfer  rows."  In  all  case3,  the  free  energy  parameters 
0-e  the  natural  logs  of  the  several  pK's  in  the  mole 
fraction  scale  which  give  the  best  fit  to  the  hemoglobin 
titration  curve  and  to  the  amount  of  carbaraino  formed 
in  the  oxygenated  and  reduced  cases. 

It  now  remains  to  show  that  the  results  of  the  con¬ 
ceptual  compartments  and  transfer  rows  in  the  model  are 
simply  mass  action  equations  similar  to  Eqs,  (15)  and 
(17);  i.e.,  that  H+  ion  and  CC^  reactions  with  hemoglobin 
have  the  hypothesized  effect  upon  the  oxygenation 
reactions . 

Write  the  mass  action  equations  from  the  model 
treating  the  transfer  coefficients  as  a  species  (see 
Table  XI).  Consider  the  first-stage  oxygenation  of  Hb^, 
the  row  "Hb402"  in  Table  XI: 


(Hb402)(RHeme)3  (OHeme) 


-15.7665 

e 


(Hb4)(02) 


(20) 
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where  the  symbol  RHeme  stands  for  a  reduced  heme  of  which 
Hb^C^  contains  three,  and  QHeme  represents  an  oxygenated 
heme.  The  concentrations  are  in  mole  fractions  and  the 
exponent  is  log  base  e  of  the  equilibrium  constant,  i.e., 
the  intrinsic  constant  Kj.  Equation  (20)  is,  then,  the 
first-stage  oxygenation  with  extra  terms  which  transform 
the  apparent  oxygenation  (as  measured  by  the  actual  con¬ 
centrations)  into  the  intrinsic  equilibrium.  Equation 
(20)  shows  this  more  clearly  if  written 


(HbA°2) 
(m>  J(o2) 


1  -15.7665 

- 3 -  e 

(  RHeme )J(0Heme) 


(21) 


where  the  concentrations  on  the  left  are  measured. 

Investigating  the  nature  of  the  transformation  on 
the  right  side  of  Eq.  '21)  via  the  mass  action  equations 
for  the  chemical  equations  in  the  conceptual  compartments 
(Table  XI)  yields,  comparing  Eqs.  (15): 


(RHeme)  (H*)  - 22. 4397 

(«2mn 


(22) 


where 
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(RHeme)  -  (HHb+  )  (the  zwitterion)  , 


K_  <«■'>(«*> 

Kl  (RHeme) 


-21.0582 
e  , 


K_  .  (Hb‘ )  (H*) 2 
^IC  (RHeme) (C02) 


-31.4723 


and  a  similar  set  from  the  oxyheme  compartment.  These 
exponents  are  in  natural  logs  and  mole  fraction  scale. 

Considering  that  the  mole  tractions  of  all  species 
in  a  conceptual  compartment  must  sum  to  one, 


(RHeme) 


1 


kri 

<H+) 


+  K 


01 


KRC(CQ2> 

(H*)2 


1 


and,  similarly, 


(OH erne) 


s& h 

“oz 


+  1  + 


01 

(H+) 


*01 


*oc(C(V 

(H+)2 


(23) 


In  addition,  the  condition  that  one  hemoglobin  yields 
four  conceptual  hemes  with  no  change  in  free  energy  gives 


Hb4  -  4(RHeme)  .  (24) 

Substituting  Eqs.  (23)  and  (24)  into  Eq.  (21)  immediately 
results  in  Eq.  (17).  A  similar  development  yields  Eq. 
(19).  The  ancillary  reactions  in  the  model  thus  affect 
the  oxygenation  reactions  in  the  hypothesized  way,  i.e., 
Eq.  (19). 
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Chapter  III 

VALIDATION  OF  THE  MODEL 

1.  INTRODUCTION 

The  addition  of  the  serum  albumin  buffering  and 
hemoglobin  reactions  permits  further  testing  of  the  mathe¬ 
matical  model  begun  in  Chap.  I.  These  tests  will  subject 
the  improved  model  to  varied  chemical  stresses  and  note 
in  detail  the  qualitative  and  quantitative  responses. 
Generally,  these  computer  validation  experiments  will  be 
similar  to  laboratory  experiments,  so  that  results  in  the 
two  facilities  may  be  compared.  Emphasis  is  therefore 
placed  on  reproducing  laboratory  conditions  as  nearly 
as  possible,  so  that  the  experiments  are  comparable. 

The  subsequent  discussions  of  particular  experiments  will 
further  illustrate  this  point. 
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2.  THE  STANDARD  STATE 

The  standard  reference  state  for  the  blood  simulation 
continues  to  be  the  arterial  blood  of  a  resting,  adult 
male,  as  in  Tables  1  and  II  of  Chapter  I,  where  various 
stages  of  an  increasingly  complex  model  approximated  these 
reference  data.  The  incorporation  of  serum  protein  buf¬ 
fering  and  hemoglobin  chemistry--the  final  stages  in  this 
series  of  experiments--requires  another  comparison  of  these 
standard  data  with  those  from  the  present  model. 

Tables  XII  and  XIII  list  the  input  for  this  model, 
called  "Standard  Model  BFFR-l."  In  comparison  to  Model  G, 
Chap.  I,  it  contains  several  important  additions.  The 
species  protein  (PROTN)  in  the  plasma  compartment  now  has 
subcategories  associated  with  the  conceptual  compartments 
Site  1,  Site  2,  Site  3,  and  Site  4,  where  the  carboxyl, 
imidazole,  c-araino,  and  phenolic  groups .  respectively, 
are  ionized.  The  miscellaneous  protein  (sometimes  called 
"X  ")  in  the  red  cells  also  is  ionized  at  an  average  pK 
of  6.79.  The  hemoglobin  is  oxygenated  in  the  four  Adair 
stages  and  undergoes  Ionization  and  carbamino  reactions 
according  to  the  data  developed  in  Chap.  II.  (Table  XIII 
also  shows  the  subcompartraents  of  the  red  cell.) 
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Table  XIII 

EXPECTED  OUTPUT  SPECIES 
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Table  XII  gives  the  total  of  species  in  the  liquid 

* 

phase  in  moles-per-liter-of-blood.  The  species  H  and  OH 
sre  combined  in  the  model  according  to  (H^)(OH 
to  give  water.  The  total  moles  of  miscellaneous  macro- 
molecular  species  and  protein  for  plasma  and  red  cells 
come  from  Dccumenta  Gelgy  [29].  Table  XII  also  shows  the 
gas  phases  for  equilibration. 

The  free  energy  p&rameterr  of  Table  XIII  resemble  those 
of  the  previous  models  (Chap.  I).  Each  additional  chemical 
reaction,  however,  also  requires  a  parameter.  For  example, 
consider  the  carboxyl  ionisation  of  serum  protein  given 
in  Site  1,  Table  XIII.  The  conceptual  compartment  Site  1 
shows  two  reactions.  The  symbol  RCOO  represents  a  carboxyl 
site  or.  a  serum  protein  *nolecule.  This  site  has  two  pos¬ 
sible  expected  states  related  by 


(RCCO~)(H~f)  . 
(RCQOH) 


U) 


where  <n  K  *  -12.40503  with  the  concentrations  in  mole 
e 

fractions,  or  -log^X  <•  *  3.64,  a  "best  fit"  average 

value  from  the  titration  data  for  carboxyl  ionisation  sites. 

Similarly,  the  remaining  ionisation  parameters  are 
derived  on  a  mole  fraction  scale  from  the  best  laboratory 
estimates  as  discussed  in  the  previous  chapter. 
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The  oxygenation  equilibrium  constants  have  the 
(Roughton)  ratios  derived  in  the  previous  chapter.  Table 
XIII  expresses  them  in  the  mole  fraction  scale  for  all 
constituents  of  the  mass  action  equations.  While  the 
ratios  of  the  oxygenation  constants  are  fixed,  the 
absolute  values  are  adjusted  empirically  to  gi^e  the 
correct  oxygenation  at.  p50,  the  oxygen  pressure  at  half¬ 
saturation.  This  empirical  fitting  is  the  best  procedure 
presently  available. 

Table  XIV  prints  the  computer  output  for  the  standard 
arterial  state.  All  species  appear  in  moles,  mole  frac¬ 
tions,  and  moles-per-liter  of  water  in  each  compartment. 

Similarly,  Table  XV  is  a  computer  solution  for  the 
distribution  of  species  in  standard  venous  blood.  The 
differences  between  arterial  and  venous  result  entirely 
from  a  change  in  gas  pressures.  For  arterial  blood,  pC^  * 
100  mm,  pCO^  *  40  mm;  and  for  venous,  pO^  *  40  mm,  pCO^  » 

46  mm.  Table  XVI  summarises  these  data  for  some  of  the 
important  species,  and  compares  them  with  similar  data 
from  the  reference  literature. 

Table  XVI  has  fsv  important  discrepencies .  The  mathe~ 
matical  model  uses  a  slightly  higher  (5  percent)  solubility 


r 
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STANDARD  ARTERIAL  BLOOD 
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-0. 
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constant  for  oxygen  than  did  the  calculations  of  the 
standard  reference. f  This  slightly  increased  solubility 
causes  a  small  rise  in  the  saturation  of  hemoglobin,  but 
the  absolute  level  of  the  oxygenation  constants  for 
hemoglobin  are  adjusted  empirically  to  fit  the  saturation 
curve.  The  carbamino  CC^  is  also  more  abundant,  in  ac¬ 
cordance  with  what  are  believed  to  be  better  laboratory 
* 

data.  The  hematocrit  is  computed  using  volume  factors 
(cc  of  water  per  cc  of  plasma  cr  red  cells)  of  94  percent 
for  plasma  and  72  percent  for  red  cells. 

In  addition  to  comparing  the  data  of  Table  XVI  with 
the  literature,  many  other  criteria  exist  for  evaluating 
a  model  of  blood  chemistry.  From  Table  XIV,  one  can 
compute,  for  example,  the  ratio  of  Cl  ^  ion  mole- fractions 
between  plasma  and  red  cell, 


iCl  ) 


PL 


(Cl*) 


1.94 


1.4s) 


RC 


(2) 


while  Bernstein  30'  gives  1.44  for  this  ratio.  Since 


The  present  data  are  from  F..I.W.  Houghton.  F.R.S., 
Cambridge  University,  private  correspondence.  19hS. 
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m 

Cl"  ions  are  assumed  not  actively  pumped,  the  Cl  gradient 
is  a  measure  of  the  Gibbs-Donnan  gradient  and  of  the 
specif Lc  ion  potential  across  the  cell  membrane.  In  this 
connection,  as  in  Chap.  I, 


(Cr)FL  t  (hco-)pl  (<)rc  (h+)rc  (oh-)pl 
(c1‘)rc  <hco^rc  (K>pl  <h+>pl  <oh'>rc 


for  the  ratio  for  all  single-valent  ions  having  the  iden¬ 
tical  apparent  level  of  the  active  pump  mechanism.  This 
study  assigns  a  zero-level  active  pump  to  species  satis¬ 
fying  Eq.  (3),  and  names  these  ions  the  "free"  ions  in 
contrast  to  ions  associated  with  a  non-zero- level  active 
pump.  The  simpler  models  of  Chap.  I  could  also  show  this 
ratio  equal  to  the  square  root  of  the  sulfate  ratio,  but 
it  does  not  hold  here,  as  will  be  demonstrated  below.  The 
constant  before  the  bicarbonate  ratio  is,  here,  1.15.  The 
HCO^  ratio  is  not  the  same  as  the  Cl”,  owing  to  the 
increased  solubility  of  red  cell  CO^  over  plasma  and  a 
slight  difference  in  apparent  pK  for  the  red  cell  bicarbo¬ 
nate  reactions.  This  Memorandum  thus  postulates  a  zero- 
level  active  pump  for  the  bicarbonate  ion  and  accounts 
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for  the  anomalous  bicarbonate  ratio  by  differences  in 
intracellular  ionic  activity--although  these  two  phenomena 
may  be  difficult  to  distinguish. 

One  of  the  major  discrepancies  in  this  model  of  human 
bloc  1  arises  from  an  inconsistency  in  the  variable  bio¬ 
logical  data.  Whereas  the  literature  [29, 30'1  gives  the 
Cl  ion  concentration  ratio  as  1.44,  it  usually  gives  the 
H  ion  concentration  ratio  at  around  1.57.  This  H  ton 
ratio  is  computed  from  the  difference  in  pH  between  plasma 
(pH  7.39)  and  red  cells  (pH  7.19),  a  difference  of  0.2  pH 
[7,311.  If  neither  of  these  ions  is  actively  pumped, 
the  cause  of  the  difference  in  concentration  ratios,  if 
it  exists,  is  not  clear.  Of  course,  a  change  in  the 
activity  of  H+  ion  or  of  Cl  ion  is  possible,  e.g.,  by 
interaction  with  the  proteins. 

In  Table  XIV,  the  model  has  the  "correct"  Cl  ion 
concentration  ratio  [30].  In  Table  XVII,  the  model  has 
the  "correct"  H+  ion  ratio  ’7\  In  each  table,  the  other 
ion  of  the  twp  has  an  e. ror  in  the  quoted  ratio.  The 
change  in  gradient  of  the  free  ions  was  obtained,  in 
this  case,  by  increasing  or  decreasing  the  average  nega¬ 
tive  charge  on  the  serum  proteins  by  4  percent,  the  lower 
average  negative  charge  giving  the  higher  chloride  ratio. 


mu****  WiW  •»»■-»  liN-OHAW  *•* 
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This  small  inconsistency  has  not  been  resolved.  The  model 
with  the  correct  chloride  ion  was  chosen  for  farther  com¬ 
putation,  since  it  also  has  the  "correct"  hematocrit  r7]. 

Foi  the  free  ions,  Eq.  (3)  holds  and  permits  the 
following  definition  of  specific  free  ion  potential: 


RT  *  n 


<C1~>PL 

<C1‘>RC 


RT (0 . 3716) 


z  F  E  , 
C  l‘ 


w 


where : 

R  *  Gas  constant  *  1.987  cal  deg  *  mole 

T  -  Absolute  temperature  *  310°, 

z  **  Valence  of  ion  *  -1, 

Cl' 

*  - 1 
F  -  One  baraday  »  23,ObO  cal  volt  mole 

E  »  Specific  ion  potential; 


or 


blS. 9(0. 3 Tib) 
(  1) (23.060) 


-  -9.9  millivolts 


(S) 


for  the  rod  ceil  free  specific  ion  potential. 
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The  specific  ion  potential  as  measured  by  Cl  ion 
gradient  differs  l.rom  that  which  could  be  computed  simi¬ 
larly  from  Na+  or  K+  gradients,  because  of  the  active 
cation  pumps.  The  direction  of  the  observed  Na+  gradient, 
and  hence  the  mode?,  gradient,  is  opposite  to  that  given 
by  the  Cl  specific  ion  potential.  specific  ion 

potential  would,  in  fact,  concentrate  the  positive  ions 
in  the  interior  of  the  cell  (see  Chap.  I),  just  as  the 
negative  Cl  ion  is  concentrated  externally.  Therefore, 
the  Na  pump  must  overcome  the  (Dcnnan)  electrostatic 
gradient,  as  well  as  raise  Na+  ions  up  to  the  concentra¬ 
tion  r^tio  observed.  Analytically, 


(Na+)pL  (Cl') 

■  ■  X  '  -  ■ 

(Na\C  <C1_)RC 


3.69 


1.45 


8.25 


2.113 

e 


(6) 


where  the  exponent  on  the  right  is  exactly  the  free  energy 
parameter,  in  Table  XIII,  for  the  Na+  pu-p.  Taking  the 
log  of  both  sides. 


i  n 


<x,\l 

<s*\c 


(Cl  ) 


-  t  n 


PL 


^  >RC 


2.  ii 


<  1.739  . 


(7) 
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Or,  if  we  write 


AF' 


RT  *n 


(8) 


The  interpretation  of  Eq.  (8)  in  this  context  is  actually 
not  .s  clear  as  Eq.  (9)  would  suggest.  Thus  £F'  replaces 
the  usual  AF,  because  the  ideal  equation  (9)  gives  only 
the  theoretical  work  required  to  move  an  infinitesimal 
increment  of  Na+  ions  (plus  the  concomitant  negative  ions) 
from  one  concentration  to  the  other  while  maintaining  the 
system  in  the  present  state;  but  it  Ignores  any  energy 
losses  in  the  membrane,  the  concentration  effects  in  each 
phase,  or  any  mechanism  which  maintains  the  gradient. 
Particularly,  Eq .  (9)  does  not  measure  the  work  required 
to  obtain  the  present  iistr ihut ion  from  classical  equi- 


libr  i  urn . 
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•f 

The  K  ion  is  concentrated  In  the  cell  interior,  so 
the  Gibbs-Donnan  gradient  aids  its  active  pump.  There¬ 
fore  , 


(K+> 


■ir. 


PL 


"4* 

(K+) 


-3.0220  -  f  n 


(Cl  ) 


PL 


PC 


(Cl  ) 


-3.3943  , 


(10) 


RC 


where  -3.0220,  from  Table  XIII,  is  the  free  energy  param¬ 
eter  representing  the  work  per  mole  to  maintain  the  ob¬ 
served  potassium  gradient. 

Computing  the  ;'potential"  of  the  K+  ion  from  Eq .  (4) 

gives 


(K+) 


ET  ’  n 


PL 


00 


RC 


613.9{-3.3943) 

(+1)23.060 


*  -9.0?  m lilt vo its  ; 


and  for  Na 


£ 
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615.9(1,798) 
Na+  ”  (+1)23,060 


-  +4.64  millivolts  . 


However,  neither  of  these  last  two  "potentials"  are 
significant  for  the  activity  of  the  cell,  since  the  cation 
mobilities  are  so  low  in  relation  to  those  of  the  anions. 
For  example,  the  "membrane  pocential"  ot  the  cell,  given 
by  Ref.  32, 


E 

m 


RT  ,  VK\c  *  PN-a<Na\c  +  PC1<C1'>PL 
F*  '  Pk(K+)  +  ?MJNa+)DI  +  Pr1  (Cl*), 


Na 


PL 


Cl 


RC 


(ID 


where  the  are  permeability  coefficients,  reduces  to 
just  Eq.  (4)  since  PNa  <  PR  «  P^. 

Tables  XIV  and  XV  show  a  somewhat  misleading  distribu- 
tion  of  the  species  SO^  and  HPO^,  as  well  as  Ca  and  Mg  * 
These  species  represent  the  so-called  acid-soluble  phos¬ 
phates,  the  non-protein  sulphate,  and  the  double-valent 
cations  obtained  by  precipitation.  While  the  distributions 
of  these  species  In  the  model  correspond  to  results  quoted 
in  the  literature  (cf.  Table  XVI),  the  "distribution"  as 
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such  is  not  well-defined.  These  species  are  much  more 
likely  to  combine  in  various  forms,  as  well  as  bind  to 
the  protein,  than  are  the  singly-charged  species  of  either 
sign.  But  these  reactions  have  not  yet  been  incorporated 
into  the  present  mathematical  model,  which  imposes  theo¬ 
retical  concentration  gradients  on  these  ions  to  create 
in  the  two  compartments  the  moles  of  each  species  observed 
in  chemical  analysis. 

While  Tables  XIV  and  XV  give  a  detailed  distribution 
of  species  for  arterial  aid  venous  blood,  Table  XVIII 
graphically  presents  the  differences  between  arterial  and 
venous  as  percent  change  from  the  arterial  standard  for 
both  the  plasma  and  red  cells.  These  data  present  several 
interesting  facts  requiring,  for  explanation,  complex 
intercompartmental  relationships.  For  example,  in  changing 
to  venous  blood  from  arterial,  the  hydrogen  ion  concentra¬ 
tion  increases  proportionately  more  in  the  plasma  than  in 

4*  4  44  44 

the  red  ceils,  whereas  Na  ,  K  ,  Ca  ,  and  Mg  move  from 
red  cells  to  plasma,  and  Cl  moves  from  plasma  to  red 
cells.  Since  in  the  model  the  gas  pressure  change  alone 
Initiates  this  pattern  of  movement,  i.e.,  pO^  changing 
fro®  100  mm  to  40  and  pCO^  changing  from  40  mm  to  46,  it 
requires  explanation. 
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Table  XVIII 
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Table  XVIII  lists  concentrations  of  species  computed 
in  millimoles  per  liter  of  l^O,  although  the  plot  is  in 
percent  change  in  mole  fractions.  In  either  case,  the 
change  in  the  water  concentration  appears  as  zero  because, 
even  on  the  mole  fraction  scale,  the  molality  of  the  solu¬ 
tion  (and  consequently  the  mole  fraction  of  1^0)  is 
practically  constant  and  is  the  same  in  each  compartment. 
Water  moves,  of  course,  under  the  influence  of  various 
stresses  applied  to  the  blood;  when  ions  cross  the 
membrane,  a  proportionate  amount  of  water  goes  along  to 
maintain  the  equal  osmolality.  The  changing  concentra¬ 
tion  of  fixed  protein  in  both  compartments  shows  that 
water  moves  from  plasma  to  red  cells  in  the  change  from 
arterial  to  venous;  the  concentration  of  plasma  protein 
increases  and  that  of  red  cell  protein  decreases.  Tables 
XIV  and  XV,  where  the  water  movement  is  also  evident, 
demonstrate  this  more  clearly. 

Blood  plasma  Is  less  buffered  than  the  red  cell 
milieu  (evident  also  from  the  Astrup  diagram  studies  of 
separated  fractions  33’).  When  the  pCO,,  rises,  the  pH 
wf  the  plasma  shifts  proportionately  more  than  that  of 
red  ceils.  Thus,  as  the  CO^  accumulates  in  the  blood  and 
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(via  carbonic  anhydrase  Ln  the  red  cell)  forms  HCO^  and 
H+,  the  hydrogen  ion  is  buffered--principally  by  protein-- 
while  the  bicarbonate  ion  redistributes  according  to  the 
Gibbs-Donnan  relations  and  in  indirect  response  to  the 
cation  pumps. 

In  the  system  model,  the  bicarbonate  ion  is  just 
another  permeable,  single-valent,  negative  ion,  though 
subject  to  intraphase  reactions.  Hence,  it  finally  ac¬ 
quires  the  gradient  of  Eq .  (3),  i.e.,  before  and  after 
the  CO^  accumulation.  In  this  view,  the  apparent  Cl 
shift  (Hamburger  shift)  is  merely  a  redistribution  of  the 
total  negative  ion  pool  in  response  to  an  additional  in¬ 
crement  of  permeable  negative  ion  whose  corresponding 
positive  ion  is  buffered.  Most  of  the  additional  HCO^ 
ion,  however,  is  found  in  the  red  cell  because  the 
increment  in  acidic  pH  shift  is  smaller  than  that  of 
plasma.  By  the  Henderson-Hasselbach  equation,  the  product 
(H*)(HCO^)  is  directly  proportional  to  pCO^;  but  with  a 
positive  increment  in  pCO^,  if  the  (H4)  does  not  increase, 
then  the  (HCO^)  must.  Thus,  as  Tabl»ss  XIV  and  XV  demon¬ 
strate,  a  large  proportion  of  the  produced  HO,  stays  in 
Che  red  cell  where  the  H+  ion  is  highly  buffered. 
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Actually,  these  statements  depend  more  completely 
upon  the  intercompartmental  mechanisms  than  has  been 
suggested.  Fixed  proteins  in  the  red  cell  and  plasma  plus 
the  active  cation  pumps  are  the  principal  determinants  of 
electrolyte  and  water  distribution,  because  of  the  inter¬ 
play  between  the  fixed  charge  on  the  protein  and  the 
active  pumps.  As  the  pCC^  changes,  protein  buffers  the 
H+  ion,  altering  the  net  fixed  charge.  Now  a  propor¬ 
tionate  change  in  the  charge  on  the  plasma  versus  the 
red  cell  protein  exists  for  which  the  Dorman  gradient 
would  remain  constant,  but  such  an  exact  proportional 
change  would  be  coincidental.  Generally,  in  the  physio¬ 
logical  range,  the  hemoglobin  buffers  better  than  serum 
albumin;  the  charge  increment  thus  differs  for  the  two 
proteins,  and  ions  must  move.  As  ions  move  across  the 
membrane,  taking  along  osmotic  water,  the  cation  gradients 
also  shift  (the  Na*  gradient  increasing  «md  K+  decreasing) 
since  the  apparent  activity  of  the  cation  pumps  is  a 
function  ot  the  Ion  reservoirs  from  and  into  which  the> 
transport  ions. 

In  addition,  the  hemoglobin  oxygenation  changes 
slightly  with  pH  and  pCO,  (the  "Bohr  effect"),  causing 
a  slight  shift  of  all  other  species.  Table  XVIII  also 
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shows  the  changes  in  the  oxygenated  species  of  hemoglobin 
in  venous  b  3od  relative  to  arterial  blood,  although  these 
changes  exceed  the  graph  boundaries.  The  fully  oxygenated 
species  Hb^Oy  decreases  considerably  while  the  remaining 
species,  insignificant  at  high  saturation,  increase  greatly 
at  low  p0?. 

THE  MATRIX  OF  PARTIAL  DERIVATIVES 

Blood  in  a  normal  unchanging  steady  state  is  probably 
very  rare.  Blood  is  a  viable  system,  so  that  even  when 
held  in  a  bottle  it  changes  chemical  composition  as  it 
metabolizes.  In  the  body,  it  also  has  a  constantly  chang¬ 
ing  chemical  composition,  at  least  cyclic  with  the  cir¬ 
culation  time.  Blood  composition  reflects  even  small 
inputs  or  losses  around  the  vascular  tree,  since  every 
species  in  the  blood  is--at  least  via  some-  series  of  mass 
action  equations--related  to  every  other. 

We  may  nevertheless  assume  instantaneous  steady 
states  at  various  and  particular  points  in  the  circula¬ 
tion  and  compute  the  chemical  composition  at  these  points. 
In  addition,  given  a  particular  steady-state  compos i t ion . 
and  given  an  increment  in  any  one  of  the  input  components, 
the  expected  response  of  each  output  species  is  comput¬ 
able.  If  b-,  i=l,...,m  are  the  input  moles  of  species, 
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and  j»l,...,n  are  the  output  moles,  then  we  compute 
the  partial  derivative,  bx^/db^,  for  each  i  and  all  j's. 

Tables  XIXA  and  B  compile  these  partial  derivatives. 
Tables  XIXA  show  the  expected  change  in  the  output  species 
(left  ;olumn)  per  unit  change  in  the  input  species  (column 
headings).  Similarly,  Table  XIXB  gives  the  expected 
change  in  total  moles  in  each  compartment.  Generally 
speaking,  small  input  increments  cause  additive  responses, 
i.e.,  the  response  to  two  small  simultaneous  inputs  equals 
the  sum  of  the  two  responses  when  the  inputs  are  added 
separately  to  the  same  reference  state.  Consequently, 
this  matrix  may  be  used  to  predict  the  systematic  re¬ 
sponse  to  small  mixed  inputs--e . g . ,  the  addition  of  NaCl 


in  water  solution, 


As  an  example  of  the  matrix's  use,  first  consider 
the  addition  of  1  meq  of  Na+  Ion  to  a  liter  of  arterial 
blood.  Reading  under  the  Na*  column,  9.13  *  10  ^  meq  ot 
the  added  Na  will  appear  in  plasma  and  the  remainder, 

_  n 

8  65  *  10  meq,  will  appear  in  red  cells.  The  addition 

- 1  -  -  2 
also  causes  5.99  »•  10  meq  of  Cl  ion  and  2.06  *  10  mM 

of  water  to  move  from  the  red  cells  to  the  plasms.  Every 

other  species  will  also  change.  H*  loti  decreases  in  both 

compartments,  but  In  the  red  cell  compartment  the  OH*  ion 
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Table  XIXA 

OUTPUT  SPECIES  INCREMENTS  PER  UNIT  INCREASE 
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Table  XIXA--conf,inued 
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Table  XIXB 

INCREMENTS  IN  TOTAL  MOLES  IN  EACH  COMPARTMENT 
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decreases  considerably  more.  In  response  to  Na+  alone, 
then,  the  red  cells  would  go  very  acidic  and  the  plasma 
alkaline.  Also,  partly  in  response  to  the  pH  change, 
hemoglobin  oxygenation  decreases. 

This  experiment  cannot  be  performed  in  the  laboratory 

since  a  negative  ion  must  also  be  added.  If  1  meq  of  Cl 

is  a:  ed  simultaneously,  however,  the  total  response  is 

the  sum  of  the  two  responses  shown  in  the  columns  headed 

Na+  and  Cl  .  Thus,  there  is  still  a  net  Na+  transfer  to 

plasma;  and  although  almost  two-thirds  of  the  added  Cl 

remains  in  the  red  cell-- including  the  Cl  transferred 

in  the  Na+  column  (5.99  x  10  ^  meq)  and  the  distribution 

-  - 1 

of  the  added  Cl  (6.44  x  10  meq  in  plasma) -- there  is  a 
net  increase  of  9.57  x  10  ^  meq  of  Cl  in  the  plasma 
compartment  as  well  as  a  net  transfer  of  water  to  plasma. 

These  predicted  effects  will  be  proved  when  dry 
NaCl  is  added  as  an  experimental  stress.  Other  details 
of  the  response  to  NaCl  can  be  predicted  from  the  matrix 
as  well  as  the  response  of  the  system  to  other  small 
stresses.  For  large  stresses,  such  calculations  must  be 

*  In  the  doctoral  thesis  of  Eugene  Kagnler .  M.O., 

Temple  University  Medical  Cent.r,  Department  of  Physiology, 
Philadelphia  (unpublished). 
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validated,  since  the  partial  derivatives  are  valid  for 
only  that  system  state  for  which  they  were  calculated. * 

OXYGEN  DISSOCIATION  CURVES 

Any  simulation  of  the  respiratory  chemistry  of  blood 
should  reproduce  the  standard  oxygen  dissociation  curves-- 
the  oxygen-loading  and  -unloading  functions  of  normal 
blood--when  the  pressure  of  oxygen  to  which  it  is  sub¬ 
jected  varies  from,  say,  1  to  100  mm.  This  test  is 
fundamental . 

Since  1942,  the  data  of  Dill  [161  on  oxygen  dis¬ 
sociation,  from  the  Mayo  Aero-Medical  Unit  have  been  the 
best  available  for  human  blooa.  However,  improved  tech¬ 
niques  (the  pH  and  oxygen  electrodes)  and  more  csrefu1 
attention  to  detail  (especially  at  very  low  and  ver\  high 

* 

hemoglobin  saturation)  may  prod  _c  an  improved  curve  ’ 34 1 . * 
The  new  data  may  modify  the  presently  best  available 
saturation  curve  by  as  much  as  2  percent  of  its  value  in 
some  regions  of  p0,?.  This  amount  is  negligible  for 
clinical  applications  where  calculations  are  seldom  within 

♦ 

See  the  discussion  of  Model  F  in  Chap.  I,  Sec.  3, 

p.  31. 

i 

F.J.W.  Roughton,  F.R.S.,  private  communication, 

1965. 
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5  percent;  but  in  biochemistry,  a  2  percent  saturation 
error  in  the  curve-fitting  procedure  may  modify  the 
oxygen  association  constants  by  100  percent  or  more,  and 
lead  to  an  entirely  different  hypothetical  microchemistry, 
particularly  in  the  associated  reactions  of  hemoglobin. 

The  present  model,  as  Ref.  2,  uses  the  Dill  data  for 
validation.  The  solid  lines  of  Fig.  2  are  the  Dill  curves 
for  human  blood;  the  plotted  points  are  various  experiments 
with  the  model.  The  points  marked  by  circles  result  from 
a  series  of  computer  experiments  (described  in  detail 
below)  wherein  the  pC0?  was  determined  empirically  to  give 
the  stated  pH  of  the  Dill  curves  at  100  mm  of  0^.  These 
points  reasonably  fit  the  laboratory  data--with  less  than 
1  percent  rrcr.  This  order  of  error  mav  not  even  be 
detectable  clinically,  but  may  be  crucial  for  a  research 
hypothes is . 

The  ’’normal"  curve  (pC0?  »  40  mm,  plasma  pH  *  7.4) 

4. 

of  Kig.  2  is,  of  course,  generated  with  all  "normal" 
values  in  the  liter  of  simulated  blood,  e.g.,  normal  buffer 
base  hemoglobin,  and  gas  pressures- - a  1 1  variables  with 
standard  values  frorr  the  literature.  The  circle  data 
points  shov  the  standard  computer  results  holding  pH 
constant  over  the  range  of  p0_, .  The  laboratory  report  *16' 
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Oxygen  dissociation  curves  for  human  blood 
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does  not  Indicate  whether  the  pH  and  pCC^  were  monitored 
throughout  the  range  of  pO^  in  the  wet  experiments.  The 
maintenance  of  both  these  parameters  is  now  known  to  be 
critical  in  obtaining  an  accurate  saturation  curve. 

Dill  computed  pCO,,  for  given  pH  values  and  vary¬ 
ing  pC^,  assuming  constant  HCO^  content.  In  fact,  the 
HCO^  varies  uniformly  with  pH  which,  in  turn,  varies  with 
pt^.  Consequently,  to  maintain  a  constant  pH  over  the 
inge  of  p(>2  in  the  saturation  curve  experiment,  either 
the  PCO2  or  the  non-carbonate  base  must  gradually  change 
to  hold  pH  and  HCO^  constant.  If  pCC^  is  held  constant, 
either  NaOH  or  HC1  must  be  added  to  hold  the  pH. 

Conversely,  to  obtain  the  extreme  pH  ranges  of  7.2 

and  7.6  indicated  in  Fig.  2  by  varying  pCC^  only,  pCC>2 

values  different  from  those  quoted  by  Dill  are  recessary. 

Whereas  D^il  computed  a  pCC^  of  25.5  for  pH  7.6  and  61.3 

for  a  pH  of  7.2,  the  corresponding  pCC^  values  are  in 

fact  closer  to  20  mm  and  85  mm,  respectively.  The 

Henderson-Hasselhach  equa' 'on  also  verifies  these  latter 

figures,  when  HCO^  varies  in  the  presence  of  a  protein 

* 

buffer  having  pK  7.2.  The  crossed  points  ly<ng  on  or 
See  F.dsall  and  Wyman,  Ref  8,  p.  569. 
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near  the  corresponding  Dill  curves  were  obtained  by 
changing  the  pCC^  to  these  new  values  and  recalculating 
the  saturation  curves  at  pH  7.2  and  7.6. 

For  all  three  of  the  curves  represented  by  the 
crossed  points,  the  computer  results  are  satisfactory, 
i.e.,  practically  indistinguishable  from  the  Dill  curves. 
This  is  important  to  several  pertinent  theoretical  ques¬ 
tions.  For  example,  the  interaction  of  and  H+  with 
oxygenation  is  evidently  simulatable,  if  not  in  exact 
detail  than  in  overall  effect.  Again,  note  that  the 
oxygenation  parameters  were  determined  for  hemoglobin  in 
solution  and  not  for  whole  blood  [21.  Yet  when  the 
parameters  were  incorporated  in  the  whole  blood,  they 
gave  the  same  curve;  i.e,,  probably  no  unaccountable 
change  occurs  in  the  oxygenation  of  hemoglobin  incorporated 
into  whole  blood  compared  to  hemoglobin  solution. 

A  further  Implication  of  the  goodness  of  fit  concerns 
the  effect  of  pH  on  the  saturation  curve.  From  p09  -  100  mm 
to  p0?  -  1  mm,  the  pH  along  the  saturation  curve  drifts 
basic  by  about  0.1  pH'  unit.  (Th’ss  amount  of  drift  is  less 
than  that  predicted  from  hemoglobin  solution  experiments-- 
0.21  pH  units  14  --because  the  blood  is  more  highly  buf¬ 
fered.)  However,  the  reasonable  fit  to  the  Dill  curve 
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indicates  that  the  pH  was  not  controlled  in  the  wet 
experiments,  and  must  have  drifted  approximately  the 
same  amount. 

With  constant  non-carbonate  base  in  the  blood  and 
constant  pCC^ ,  the  pH  of  normal  blood  goes  basic  when 
the  blood  deoxygenates,  and  vice-versa  [141.  The  exact 
amount  of  this  effect  varies  both  with  the  absolute 
saturation  and  the  absolute  pH.  Beginning  at  100  mm  pC^ 
and  pH  7.4,  deoxygenation  causes  the  non-ccnstant  pH 
curve  to  deviate  from  the  constant  pH  curve,  but  by  an 
amount  different  than  for  the  7.2  or  7.6  curves  [351. 

The  oxylabne  side  reactions  of  the  protein  produce  this 
phenomenon,  related  to  the  Bohr  ef£  ct.  Certain  ionizing 
groups,  as  yet  not  precisely  specified,  change  pK  ac¬ 
cording  to  the  state  of  a  nearby  oxygen  site.  Generally 
speaking,  the  pK  of  nearby  sites  decreases  as  oxygen 
attaches,  although  for  certain  sites  it  may  rise  with 
oxygenation  [36'.  Also,  some  of  these  sites  are  evidently 
carbamino  sites,  since  the  absence  of  CO.,  reduces  the 
shift  in  pH  effect  [35]. 

To  teat  the  difference  between  a  curve  run  at 


constant  pH  and  one  in  wliich  the  pH  drifted,  HC1  or  NaOH 
was  added  as  though  a  laboratory  pHStat.  had  been  used  in 
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the  wet  experiment.  The  circle  points  plotted  in  Fig.  2 
indicate  the  results. 

While  the  error  from  the  Dill  curves  is  biochemically 
important  at  low  pO^,  Che  absolute  amounts  of  these  effects 
are  difficult  to  verify  in  the  literature,  because  of  the 
absence  of  definitely  comparable  experiments.  In  several 
published  experiments  which  study  the  pH  effects  on 
oxygenation,  principal  results  are  on  hemoglobin  solution 
[14,28,36,37].  Since  the  whole  blood  is  buffered  more 
heavily  than  hemoglobin  solution,  the  quantitative  values 
for  whole  blood  are  not  evident  from  these  experiments. 
Several  other  papers  note  that  for  whole  bl^od  the  de¬ 
crease  in  pH  shifts  the  saturation  curve  to  the  right, 
i.e.,  decreased  oxygen  saturation  at  the  s...  j  p0?  'V3S 
expected);  but  generally  this  shift  in  pH  is  accomplished 
by  changing  pCO,,  [38,39;.  In  short,  a  careful  test  of 
this  computer  experiment  would  be  interesting. 

TV! K  !,ASTRUPM  EXPERIMENTS 

a 

Figure  3  reproduces  an  Astrup-type  chart  for  clinical 
determinations  of  the  state  of  whole  blood.  The  straight 

♦ 

Copyr  ight  *  RADIOMrTX  £R  ,  Copenhagen  . 
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Fig. 3— Astrup- type  experiments 
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line  A  (from  the  literature  '4CT)  represents  an  experi¬ 
ment  with  normal  adult  male  blood;  D  is  a  comparable 
experiment  computed  from  the  present  model. 

Generally  speaking,  the  slope  of  the  lines  drawn  on 
the  Astrup  chart  from  blood  samples  is  proportional  to 
the  buffer  power  of  the  blood  in  respiratory  acidosis. 

The  steeper  the  line,  the  better  the  blood  will  buffer 
or  increase  in  pCO^  (i.e.,  H+  ion).  Thus,  the  line  A-- 
or  the  line  D  (since  they  differ  insignif icantly)--gives 
the  norr..al  buffer  slope  for  blood  with  the  standard 
amounts  of  buffering  constituents.  In  Fig.  3,  blood  at 
pH  7.38  and  at  40  mm  pC0o  contains  the  amount  of  "standard" 

a. 

bicarbonate  given  by  the  Henderson-Hasselbach  equation, 
here  23.8  meq/L  blood.  With  other  conditions  constant, 
but  changing  pC09  ,  the  buffered  blood  will  normally  follow 
line  A  (or  D)  intersecting  the  hemoglobin  content  curve 
at  16  (17)  grams  percent. 

The  buffering  constituents  are,  in  decreasing  order 
of  impor tance ,  the  carbon  dioxide  system  itself,  hemo¬ 
globin,  plasma  proteins,  organic  and  inorganic  phosphates, 
and  subsidiary  systems.  Togethe’'.  thest  several  con¬ 
stituents  represent  the  anions  capable  of  binding  hydrogen 
ion  in  the  viable  pH  range  and.  except  for  the  CO.,  system. 
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are  called  the  "buffer  base"  of  the  blood.  The  carbonate 
system  is  not  part  of  the  buffer  base  because  it  is 
(indirectly)  the  CO.,  which  is  to  be  buffered.  "Buffer 
base"  is  an  imprecisely  defined  term,  since  the  hydrogen¬ 
binding  capacity  of  the  proteinate  and  phosphate  factors 
varies  so  considerably  with  pH,  their  "buffer  power"  being 
directly  proportional  to  the  slope  of  their  titration 
curves  at  a  particular  pH.  Nevertheless,  the  term  "buffer 
base"  is  useful  for  the  narrow  range  of  pH  in  clinical 
applications.  At  normal  concentrations  of  protein  end 
at  viable  pH  and  pCO^,  the  imidazole  ionizing  groups  of 
protein  plus  the  carbamino  reactions  and  phosphate  com¬ 
plexes  normally  buffer  as  though  about  45  to  4b  meq  per 
liter  of  ionizing  sites  existed. 

Additions  ol  acids  or  bases,  as  in  metabolic  acidosis 
or  alkalosis  without  respiratory  compensat ion ,  shouio 
change  an  available  bufter  base  of  the  blood  by  just  the 
amount  of  acid  or  base  added.  Because  the  solvent  water 
ha*  a  low  ionization  constant,  the  added  excess  of  hydrogen 
or  hydroxyl  ion  can  not  stay  in  solution  in  ionized  form; 
these  ions  must  essentially  either  bind  to  a  hydrogen¬ 
binding  site  or  remove  a  bound  hydrogen- -  this  reducing  or 
increasing,,  respectively,  the  available  buffer  base. 
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When  the  buffer  base  is  non-normal,  indicating  a 
metabolic  alkalosis  or  acidosis,  the  buffering  reactions 
of  the  blood  to  increased  pCX^  are  approximately  the  same 
as  for  the  normal.  However,  the  buffer  line  on  the  Astrup 
diagram  will  shift  to  the  new  range  of  pH  corresponding 
to  the  new  value  of  the  buffer  base.  Lines  A'  and  D'  of 
Fig.  3  are,  respectively,  the  ideal  responses  of  whole 
blood  in  the  Astrup  experiment  and  the  response  of  the 
computer  model  to  changes  in  after  addition  of  10 

meq  of  HC1  per  lite*\  The  slope  of  the  buffer  lines  has 
steepened  slightly  because  the  lower  pH  is  closer  to  the 
middle  o:  the  range  of  imidazole  ionization,  so  the 
buffering  is  more  efficient.  The  lines  D'  did  not  move 
a  full  ten  units  or*  the  base  excess  curve,  although  10 
meq  of  acid  was  added. 

This  discrepancy  is  difficult  to  rationalize.  The 
literature  notes  that  the  standard  curves  are  based  on 
blood  samples  from  12  persons,  and  that  the  same  titra¬ 
tion  results  were  obtained  with  hydrochloric,  lactic,  or 
acetic  acids  [40],  Yet  at  constant  pCO-,  the  computer 
blood  buffers  HC1  better  than  real  blood;  i.e.,  for 
constant  pCC^ ,  10  meq  of  HC1  changes  the  pH  less  In  the 
model  than  the  same  amount  added  to  human  blood. 
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The  discrepancy  is  large — about  10  percent  of  the 
buffering  capacity  in  this  central  pH  range;  and  its 
direction  indicates  that  the  model  is  already  buffering 
too  well — i.e.,  the  model  does  not  need  an  additional 
buffer  system.  On  the  other  hand,  the  blood  used  to  con¬ 
struct  the  chart  was  probably  not  viable  blood;  it  was 
probably  damaged  by  centrifugation,  changes  in  temperature 
and  gas  pressures,  and  by  the  addition  of  a  fluoride 
stabilizing  agent  [40].  Some  denaturation  of  the  protein 
by  physical  forces,  plus  slight  hemolosis,  plus  fluoride 
agents'  interference  with  cell  glycolysis  could  easily 
make  the  blood  a  10-percent- less-efficient  buffer. 

Line  B  in  Fig.  3  is  the  iso-bicarbonate  line,  i.e., 
the  line  along  which  pH  +  log  pCC^  “  K  is  constant,  and 
the  HCOj  concentration  is  consequently  also  constant. 
Ed8all  and  Wyman  [8]  give  line  E  for  the  pCC^/pH  relation¬ 
ship  in  0.02  M  sodium  bicarbonate  solution;  the  HCQj 
concentration  is  nearly  constant.  They  also  remark  that, 
for  this  simple  solution  esse,  the  total  C0n,  T t  carried 
is  relatively  invariant  with  pCO^,  the  change  being  only 
from  17.1  to  22.0  ccM  per  liter,  with  e  change  in  pCO^ 


from  0.45  to  63.0  ca. 
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Blood,  however,  for  efficient  action  must  have  a 
different  property.  Relatively  large  ar><ounts  of  CO^ 
must  be  picked  up  or  lost  with  small  increments  of  pCC^ 
between  venous  and  arterial  blood.  This  is  accomplished 
through  better  buffering  of  the  blood  than  the  simple 
solution  affords.  If  the  pH  is  almost  constant,  then,  by 
the  Henderson-Hasselbach  equation,  as  pCO^  increases 
HCO^  must  increase  proportionately.  Thus,  buffering  the 
hydrogen  ion  is  essential  to  CO^  transport. 

From  an  earlier  mathematical  model  of  human  blood  [ 2 1 
not  containing  the  imidazole  reactions  in  plasma  and  red 
cell  protein,  line  C  in  Fig.  3  indicates  the  buffer  power-- 
and  hence  CO^-transporc  capacity.  Rotating  this  line  to 
the  present  increased  buffer  capacity  (line  D)  was  a  matter 
of  including  the  titration  reactions  of  plasma  protein  and 
of  hemoglobin.  Of  these  two  groups  of  proteins,  that  of 
the  red  cell  is  more  important  because  of  the  carbamino 
reactions  of  hemoglobin  with  C0? ,  oxyiabile  ionizing 
groups  and  the  large  group  of  non-oxylabile  sites. 
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